STABILIZATION OF THE LINEAR SYSTEM OF MAGNETOELASTICITY 
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Abstract. We give a necessary and sufficient condition, of geometrical type, for the uniform decay 
of energy of solutions of the linear system of magnctoelasticity in a bounded domain with smooth 
boundary. A Dirichlet-type boundary condition is assumed. When the geometrical condition is not 
fulfilled, we show polynomial decay of the energy, for smooth initial conditions. Our strategy is to use 
micro-local defect measures to show suitable observability inequalities on high-frequency solutions of the 
Lame system. 



1. Introduction 

1.1. The system of magnetoelasticity. Let fi be a bounded, simply connected domain of R 3 , with a 
smooth boundary. Let us consider the following system, modelling the displacement of a elastic solid in 
a magnetic field: 

d 2 v - [iAv - (A + ii)Vdivv - k rot/i AB = 0| 

(3d t h + rot rot/i - (3 rot (d t v A B) = > (t, y) G (0, oo) x fl 
111 div h = J 

v = 0, h.n = 0, curl h An = (t, y) € (0, oo) X dfl, 

where v — (vi, 1)2, V3) is the displacement vector of the solid, and h — (hi, h^, /13) the magnetic field. The 
system is located in a constant exterior magnetic field B = (5,0,0). We have denoted by A, V, div, 
curl respectively the Laplace operator, gradient, divergence and curl operators according to the space 
variable y, in the euclidian metric of R 3 . The positive constants k and j3 are coupling constants, and n 
is the external normal vector to the boundary of tt. The real Lame constants A and ^ are such that: 
A + 2/i > 0, fi > and A + /i ^ 0. 

The system Q has a natural time-decreasing energy: 

E {t) = l [ \d t v\ 2 + fi\Vv\ 2 + (\ + f i)\divv\ 2 + n\h\ 2 dy. 

When £1 is simply connected, G. Perla Menzala and E. Zuazua have showed that this energy tends to 
zero as time tends to infinity, which is a simple consequence, using La Salle invariance principle, of the 
non-existence of stationnary solution for |QJ . The goal of this paper is to give estimates on the speed of 
this convergence. 

The system Q may be seen as a coupling between the Lame system: 
(2) d 2 u - [iAu - (A + /it) Vdiv u = 0, 

with Dirichlet boundary conditions, which is a conservative system, and the following heat equation: 

f3d t g - Ag = 0. 

The decay of energy is produced by this strongly dissipative equation. From the point of view of v, the 
dissipation is caused by the coupling term: R(v) :— curl (dtv A B). 
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Let us first consider the uniform decay with respect to initial condition of the energy: 
(3) E{t) < f(t)E(0), f(t) — > 0, 

t— >-j-oo 

where / if independent of the initial condition. In this case it is easy to show, using the semi-group 
property of the equation (|TJ, that / maybe taken as a negative exponential function. 

In paragraph ll.2l we state, with a technical hypothesis on SI, a necessary and sufficient condition on the 
geometry of the problem for © to hold. When this condition is not fulfilled, there exist rays on Q, named 
B-resistant rays, along which the energy of some solutions of Q concentrates, and the dissipative term 
R(v) is very small. Indeed, when such a ray exists, there is a sequence of solutions of Q concentrating 
on the ray and which is in first approximation parallel to B. 

When there is no uniform stabilization we show (with the same technical property on SI than before), 
that solutions of JQ) decay with polynomial speed for smooth enough initial data (cf paragraph 11.3(1 . 
The speed of decay still depends on the geometry of SI. In this case, the possible existence of boundary 
.B-resistant rays (i.e. living only in the boundary of SI) of infinite life-length is the main obstacle to the 
decay. 

Before giving more explicit results, let us mention some earlier works on related subjects. As it was 
already stated, the convergence to for the energy of magnetoelasticity in a bounded, simply connected 
domain was shown by G. Perla Menzana and E. Zuazua in |11|. but their method does not give any 
information on the rate of convergence. By energy methods, Munoz Rivera and Racke 12 , Munoz Rivera 
and de Lima Santos |13j have shown the rate of convergence to be at least polynomial, in dimension 2 or 
3, but only for some precise types of domains. Andreou and Dassios [I] have examined the same system 
on the entire space R 3 , showing again polynomial decay for some initial conditions. 

The linear system of thermoelasticity has been more precisely understood. In this system, the Lame 
equations are coupled with a scalar heat equation. The dissipation is caused by the longitudinal part of 
the Lame equation (the curl- free part of v). In f) and [2], the authors give (under a spectral assumption) 
a necessary and sufficient condition on f2, of geometrical nature, for the uniform decay in dimension 2 or 
3. Namely, this decay is equivalent to the non-existence of rays, called "transversal polarization rays", 
carrying the transversal component of v (the divergence-free component), which resists to the dissipation. 
In the authors also prove the polynomial decay in dimension 2, under the same spectral assumption, 
which is namely that the operator associated to the equation does not admit any real eigenvalue. As 
shown in |llj . this spectral condition is always fulfilled for the system of magnetoelasticity in a bounded, 
simply-connected domain. 

The comparison of the two systems of thermo and magnetoelasticity show that thermoelasticity is 
slightly less dissipative (the coupling of the Lame system with the heat equation is weaker), and more 
difficult to describe, because of the non-trivial polarization of transversal waves. 

1.2. Uniform decay. Assume that dtt has no contact of infinite order with its tangents. Thus, the 
hamiltonian flow of the symbol of a dAlembertian 9 t 2 — c 2 A, which is defined locally in S*(M. x SI) (the 
spherical cotangent bundle of SI), maybe extended until the boundary of this bundle to a global flow, 
the generalized bicharacteristic flow, wich may be seen as a continuous flow on the spherical compressed 
cotangent bundle S£(ffi. x SI) (cf jZJ chap. 24.3]). We shall call bicharacteristic rays or just rays the 
characteristic curves of this flow. Such a curve 7 will be said parallel to B if its direction of propagation 
is always parallel to B and orthogonal to B if its direction of propagation is always orthogonal to B. We 
refer to section |3] for the exact definitions of S£(R x Vt) and of the generalized bicharacteristic flow. 

The Lame system J3J may be written as the sum of two wave equations known as the longitudinal and 
transversal wave equations, of respective speed cl '■= y/X + 2/i and ct '■= JJi (cf paragraph 13. 5(1 . The 
assumption A + (i ^ is equivalent to cx 7^ cy. 

Definition 1.1. One calls longitudinal ray (respectively transversal ray) any bicharacteristic ray for 
the operator df — c 2 A (respectively d\ — c^A). One calls B-resistant ray any continuous application 
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7 from an open interval / = (so,s n ) to S£(R X Q) such that there exists a finite number of reals 
so < si < ■■■ < s n such that: 

• on (Sj-i, Sj), j £ {1, n}, 7 is a longitudinal ray parallel to B, or a transversal ray orthogonal 
to B; 

• if j € {1, ...,n — 1}, j(sj) is an hyperbolic point for the longitudinal and transversal waves (cf 
paragraph 13. 1.7f) and one of the following assertions is true: 

— (L — ► T) case: 7 is a longitudinal ray on ]sj_i, Sj[, and a transversal ray on ]sj, Sj+i[; 

— (T — ► L) case: 7 is a transversal ray on ]sj_i, Sj[, and a longitudinal ray on ]sj, Sj+i[. 
(cf figure ^) 

Near Sj, 1 < j < n — 1, the continuity imposed by the definition of 7 gives a condition on the angles 
of incidence and refraction. In the case (L — > T), if we denote by aj, the angle between the longitudinal 
incoming ray and the tangent to dil in the plane of incidence, and by j3r the angle between the transversal 
outcoming ray and this tangent (cf figure ^ c), we have: 

tana^ = — , tanpy = — , 
cl c t 

(which implies at + @t — k/2). In the case (T — + L), and with similar notations, we have: 

— , tan/3 L = — . 
ct c h 

Remark 1.2. The B-resistant rays of figure^are all planar, but this is not a general property. 
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a) A transversal 
B-resistant ray 




b) A longitudinal 
.B-resistant ray 
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c) Transfer [L -> T) 




d) A boundary 
B-resistant ray 



B 



Figure 1 . Examples of B-resistant rays 
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Theorem 1. Let Q, be a bounded, simply connected domain o/R 3 , with a smooth boundary, having no 
contact of infinite order with its tangents. 

The energy of the system of magnetoelasticity in f2 decays uniformly if and only if there exists an L > 
such that every B -resistant ray on f2 is of length at most L. 

Remark 1.3. As it will be shown in the proof, the transversal rays carry the component of v which is 
orthogonal to the direction of propagation, and the longitudinal rays the component of v which is parallel 
to this direction. A _B-rcsistant ray, whose direction of propagation is orthogonal to B in the transversal 
case and parallel to B in the longitudinal case carries essentially the component of v which is parallel to 
B, thus cancelling the dissipative term: 

R(v) := curl (d t v A B). 
From this point of view, the theorem ^ is very natural. 

Remark 1.4. It is essential to assume cl ^ ct- Otherwise, the first equation in would be a wave 
equation with wave speed q = cy. Every solution of JJ) such that: 

«lt=o-l~B, d t vit=o-LB, hit=a = 

would be of constant energy. 

Remark 1.5. If £1 is not simply connected, there exists a finite dimensionnal space E of stationary 
solutions of (JIJ , whose components along v are null. The study of the decay to zero of the solutions may 
be replaced by the study of their convergence to the eigenfunctions corresponding to the space E (cf (lit 
chap. 5]). We won't develop this aspect here. 

Remark 1.6. The condition of uniform decay is not fulfilled in simple cases, like the one of a bowl, but 
is generic in the class of C°° open sets. 

1.3. Polynomial decay. Now we state a result of polynomial decay for initial data which are sufficiently 
smooth. The existence of a boundary B-resistant ray of infinite life-length is equivalent to the existence 
of a smooth closed curve T of 90, included in a plane V normal to B, boundary of a convex set of V, and 
such that on T, n is normal to B (cf figure^ d). On such a curve, B stays tangential to the boundary. 
Let H(r) be the minimal order of contact of dft with a tangent parallel to B. If such a curve T exists, 
and if the boundary has no contact of infinite order with its tangents, then: 

2 < K(r) < oo. 

If (v, h) is a sufficiently smooth solution of (QJ, we shall denote by E^>(t) the energy of order j of (v, h): 

E^\t) := ^-E{t). 
w dP w 

Let Xj be the subspace of X of all initial data of 0) such that E^ (0) is finite. It is exactly the domain 
of A' , where A is the linear operator of magnetoelasticity defined in section 

Theorem 2. Let n be a bounded, simply connected domain with smooth boundary having no contact of 
infinite order with its tangents. 

a) Assume there is no boundary B-resistant ray on ft of infinite life-length. Then: 

boo, yv ex u Vt>0, E(t) < : -BW(0). 

t ~\~ 1 

b) Assume on the contrary that such rays exist. LeT\, ...,T m be the support of this rays, and: 

K:= sup N(T m ). 

m=l..M 

Then: 

3C>0, VVoEX K , Vt>0, £?(t) < _^.£?W(0). 
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Remark 1.7. By an easy interpolation argument, one may deduce from b the polynomial decay of any 
solution with initial condition in X%: 

VK,eX x , Vt>0, E(t)< C E<U(p). 

Remark 1.8. Theorem El completes the works of J. E. Muhoz Rivera et M. De Lima Santos ^Sj which 
show, for some types of domains of R 3 , a decay in \jt for initial data in E^ 7 \ Note that the domains 
considered in their work (all of which have contacts of infinite order with their tangents) do not fall within 
the scope of our article. 

The remainder of the paper is organized as follows. In section [21 we reduce theorems ^ an d [21 to 
high-frequency observability inequalities on the Lame system (J2J. This is based on two arguments: the 
setting aside of low frequencies, which is a consequence of the non-existence of stationnary solution for 
the equation Q shown in 1TJ, and the decoupling, by simple calculations, of the two equations (the 
Lame system and the heat equation) which compose QJ. In section [21 we introduce micro-local defect 
measures (an object due to P. Gerard jHj and L. Tatar [^j, and in this particular setting to N. Burq and 
G. Lebeau 0), in order to study the lack of compactness of a sequence of high-frequency solutions of the 
Lame system. The main result of this section (apart from the existence of the measures) , is a propagation 
theorem which was stated and shown in [2] . In section^ we prove the observability inequality on solutions 
of the Lame system which implies theorem^ The method of proof is to introduce, in a contradiction 
argument, a sequence of high frequency solutions of J2J which contredicts this inequality, and to use 
propagation arguments on the defect measures of this sequence. Section is devoted to the necessary 
condition of theorem ^ and is inspired by [5] : defect measures are used to construct a sequence of 
solutions of ^ concentrating on a _B-resistant ray and contradicting an observability inequality. Finally, 
in sectional we prove by similar arguments than those of section 0] an observability inequality with loss 
of derivatives which implies the polynomial decay. 

The author would very like to thank his thesis advisor, Nicolas Burq for his invaluable help. 

2. Observability inequality for the Lame system 

2.1. Notations and preliminary results. In this subsection are gathered a few basic facts about 
equations QJ and (J2J, as well as some notations. The main results of section [21 are stated in the next 
subsection. 

If U is an open set of M 3 or R 4 we set: 

H S (U) := H S (U, C 3 ), L 2 {U) := L 2 (U, C 3 ). 
2.1.1. Magnetoelasticity. Consider the following spaces: 

H := {g € L 2 (fl), div g = in ft, g.n = in dft} 
Hi :={f eH 1 ^), / = 0in9ft} 
H' := {/ eHr\H 2 (n), curl/ An = in 9ft}, 

and the following norms: 

\\9\\ 2 H : = K ll5lll2( fi ), Il/ll«i(f2) : = ( A + /^)ll div /llL2(^) +Ml|V/||( i2(n)) 3. 

Let A be the unbounded operator on X := Hq{Vl) x L 2 (fl) x H, with domain D(A), defined by: 

A{V ):=\ -A e «o-K(curl/io) AB 

y — curl(ui AB) + icurlcurl/io / 

D(A) := (H 2 n Hi) x Hi x H' . 
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where Vq = (vq, Vi, ho) denotes an element of X. Equation Q may be rewritten: 

(4) d t V + AV = 0, V=(v,d t v,h). 
The following proposition is due to G. Perla Menzala and E. Zuazua jllj : 

Proposition 2.1. a) The operator A is maximal accretive. For any initial data Vq G X , there exists an 
unique weak solution V(t) = (v(t),dtv(t),h(t)) G C°([0, +oo[; X) of (@J) such that V(0) = Vq. Functions 
v and h are solutions in the distributional sense of the three first lines of system 

b) The energy E{t): 

(5) E{t)-=\\\V{t)\\ 2 x = \! \d t v\ 2 +fi\Vv\ 2 + (X + ^\dWv\ 2 + K\h\ 2 dy 
is decreasing. More precisely: 

(6) Vt > 0, E(t) - E(0) = ~ [ f \roth\ 2 dy. 

P Jo Jn 

c) If fi is simply connected: 

Wo e X, E{t) — ► 0. 

t — *+oo 

The assertions a) and b) are straightforward applications of the semi-group theory for the oprator A. 
The assertion c) is a consequence of the non-existence of stationnary solutions for the system. 

2.1.2. Lame system. Let us now consider the Lame system with Dirichlet boundary conditions: 

d 2 u - A e u = in R x fi 

(7) u ran = 

(u\- t= o,dtu\-t=o) = (u ,ui). 

Let X e be the space Hq x L 2 and C the unbounded operator on X e defined by: 



(8) C := 



-Id 
-A fi 



D(C) := H 2 n Hq 1 x L 2 



Taking (uq, u\) in the energy space X e , the equation Q may be written: 
(9) dtU + CU = 0, U(t) = (u,d t u). 

Proposition 2.2. The operator C is maximal and unitary. For any initial data Uq = (v,q,u\) G X e , the 
system has an unique weak solution U G C (M, X e ). Furthermore, the function u is a solution of 
in the distributional sense. At last, the energy: 

\{\\umHi + \\dtu\\ 2 L2 ) = \\\U\\ 2 Xc 

of this solution is constant. 

2.1.3. Two useful lemma. The two following standard lemma will be of great help in all this paper. The 
first one is due to the fact that f2 is simply connected (cf |17l Appendix I, lemma 1.6]): 

Lemma 2.3. The H 1 norm on H n H 1 is equivalent to the norm: \\u\\ := || curl ^2 . 

The second lemma is a elementary energy estimate on solutions of the non-homogeneous Lame system. 
If w(t) is a function with values in some Hilbert space, we set: (wq,wi) := (w, dtiv) |- t=0 . 
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Lemma 2.4. Let T > 0, W € C°((0,T),X e ) and F e L 2 ((0,T),X e ) such that: 

(10) d t W + £W = F, te(0,T). 

Then: 

/T \\W(t)\\l dt<c\ ||W(0)||i + f T \\F(t)\\xdt) > 



where C only depends on T. In particular, if: 

w e c 1 ((o,T) ) i 2 (o))nc ((o,T) l ff 1 (fl)), / e l 2 ((o,t) x n) 

dfw - A e w = f. 

Then: 

ll w llir 1 ((o,T)xn) - C (II^oIIk^o) + Il w illi2(^) + ||/Hi2((o,T)xn)) ■ 

Proof. To prove the first inequality, we may suppose Wo E D(C). The X e scalar product of 1)10(1 with W 
gives 

Re{d t W, W) Xe + Re(£W, W) Xe = (W, F) Xe 



\j t \\W(t)\\ 2 x c <W{t)\\xA\F{t)\\ 
J t \\W(t)\\x c < \\F(t)\\ Xc 



x, 



\\W(a)\\x. < \\Wo\\ x .+ / \\F(t)\\ Xe dt, s € (0,T) 
Jo 

\\W(s)\\ Xe <c{\\W \&.+s£\\Fm*x e dt } 



(the last line is a consequence of Cauchy-Schwarz inequality). Next, we bound, in the right member of 
the inequality the integral from to s by the same integral from to T, and we integrate with respect 
to s between and T, which yields the first part of the lemma. The second part is an easy consequence 
of it. □ 

2.2. Results. 

Proposition 2.5 (Uniform decay). Let fl be a smooth, simply connected, bounded domain o/K 3 . 

a) Assume that there exist T > and C > such that for any solution U of ffi) : 

(11) |M^i + |MH 2 < C (||curl(5 tW A B)f H - 1{(0>T)xn) + \\u \\ 2 L 2 + KH^) 

then the energy of solutions of the system of magnetoelasticity 0) decays uniformly with respect 
to initial data. 

b) Conversely, if the energy of solutions of PJ) decays uniformly, then there exist T > and C > 
such that for any solution of of finite energy: 

(12) IKHi, + Mil* < C f \\d t u A Allien). 

Jo 

Remark 2.6. The two inequalities 1)11(1 and 1(12(1 are indeed equivalent (by theorem^. 

Let U = (u,dtu) be a solution to © with initial data Uq € D(C N ). Set: 

JV 



Gr(«) :=^||cui-1(^ +1 mA£) 



2 

Lfj- 1 ((o,T)xn)- 
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Proposition 2.7 (A sufficient condition of polynomial decay). Let fl be a bounded, simply connected 
domain o/K 3 . Assume that there exist T > 0, C > and an integer N > 1 such that for every solution 
of the Lame system jnj with initial data: 

(13) U = {u ,ui)eD(C N ), 
the following inequality holds: 

(14) Ikoll 2 ^ + HHL < C (Q$(u) + \\u \\h + • 

Then there exists C > such that for every solution V of (0} with initial date Vq E D(A n ), and for all 
positive t, 

\\V{t)f x < ^HVoll^) 

Remark 2.8. results such as propositions 12.51 and 12 . 71 are fairly classical in this setting. To prove them, 
we shall avoid the usual abstract decoupling argument (see [0]) but rather use simple energy estimates 
on systems of magnetoelasticity and Lame. 

2.3. Uniform decay. We prove here the proposition 12.51 We first write a necessary and sufficient 
condition of uniform decay for solutions of a general dissipative equation. The second step of the proof 
consists in applying this condition to the system of magnetoelasticity, furthermor decoupling it in the 
system of Lame and an heat equation. 

2.3.1. Abstract framework. Let V be a maximal, accretive operator on an Hilbert space X, with dense 
domain D{V). Denote by ||..|| the norm of X, \\..\\i the natural norm of DCP 1 ) and the norm of 

its dual space, with respect to the pivot space X. Assume the embedding: 

X — > D(V)' 

is compact. For zq e X, we will denote by z(t) the solution (obtained for example by standard semi- group 
theory) of: 

(15) lti +Vz = 0, Zt =° = z ° 

By accretivity of V, the energy ^H^H 2 is time-decreasing. The following uniqueness-compactness argument 
is by now classical (cf [3]): 

Lemma 2.9. The two following assertions (i) and (ii) are equivalent: 

(i) 3C > 0, 3a > 0, Vz £ X, Vt > 0, ||z(t)|| 2 < C\\z \\ 2 e- at 
(the energy is uniformly decreasing) 

(ii) a) 3T >0,3C> 0, Vz £ X, \\z(T)\\ 2 < C (||z(0)|| 2 - ||z(T)|| 2 + ||^C0)1I— l) 
b) There is no non-zero solution of }15\) of constant energy on [0, +oo[. 

Corollary 2.10. The energy of ^j) is uniformly time- decreasing if and only if: 



3T > 0, 3C > 0, Wo E X, E(T) < C { / / |curl/i| 2 dydt 



Wo\\d{a)> 



Indeed, the non-existence of stationnary solution (the condition ii,b of lemma |2~§I has been proved in 
[111 p.356]), which shows the corollary. 

Proof of lemma Wlk It is easy to see that (i) may be replaced by: 

(f) 3T > 0, 3C > 0, Vz G X, ||z(T)|| 2 <C(||z(0)|| 2 -||z(T)|| 2 ) 

Clearly (i') implies (ii). 
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Assume (ii). For some T > 0, set: 

q T (z) := ||z(0)|| 2 - ||z(T)|| 2 , G T := {z e X, q T (z) = 0}, 

which is the kernel of a positive, bounded, quadratic form on X, thus a closed subspace of X. 

According to (ii), a), and the compactness of the embedding from X to D(P)' , G(T) is locally compact 
thus of finite dimension, for large T. By assumption b), 



n gt = {°}- 



T>0 

Consequently, dimG T being a time decreasing function of T, when T is large enough: 

(16) G T = {0} 

Let's fix such a T. The quadratic form qr is positive definite so that its square root y / gr is a pre-hilbertian 
norm on X, bounded from above by the natural norm of X. Assume (if) does not hold. Then there 
exists a sequence (z$) of elements of X such that: 

(17) l = ||z fc (T)|| 2 , lim q T (z k ) = 0. 

k — >+OG 

This implies that ||zq|| is bounded. Thus, we may extract form (zg) a subsequence, which we will again 
denote by (zg), such that: 

Zq — >■ zq £ X, weakly in X. 

k — > + oo 

Let (fT be the hermitian product given by qx- We have: 

lim ip T (z k ,z) — > q T (z), 

k — >+oc k — >+oo 

which implies, with \17j. that qr{z) = and thus, using (|16fl that z = 0. The compactness of the 
embedding of X in D{P)' yields: 

hm ||zo fc ||-i=0. 

k — >+oo 

Using a) and l|17|l we obtain the following contradictory assertion: 

1 < C (q T (z k ) + M^o II— x) = "(I) quand + oo. 

□ 

2.3.2. Proof of vrovosition [KR Assume the uniform time-decay of the energy of solutions of Then, 
by ®, there exist T > and C > such that the following estimates hold for any solution v of |J2}: 

(18) IKIIx<C f T \\cmlh(t)\\ 2 LHn) dt. 

Jo 

Let U be a solution of the Lame system with initial data Uq — (uq, ui) € D{C) and V the solution of the 
system of magnetoelasticity with initial data: 

V = (v,d t v,h)\t= = (u o ,ui,0). 

Set: W(i) := V(t) - (u(t),d t u(t),0). Then: 

d t W + AW = (0, 0, -curl (d t u A B)) 

Take the scalar product in X with W of the two side of this equality, then incgrate the real part with 
respect to time between and T. Using: 

R G {AW,W) x = j 3 \\cur\h\\l 2{n) 



(curl (d t u A B), h) L ir n \ = {d t u A B, cm\h) L2 



(O) 
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the fact that W\t=o = 0, and in the second line, the inequality (TSJ, we get: 

\\W(T)f x + f ||curl/»(t)||i, (n ) dt < C f \\d t u{t) A Bf LHn) dt 
Jo Jo 

WMh^o) + K||l 2(n) < C f \\d t u(t) A B||| 2(n) dt 

J 

This shows point b). To prove a), assume that inequality holds. Consider a solution V = (v, dtv, h) 
of |@J with initial data Vo — (vo,v\,ho), and the solution u of Lame system with initial data: 

(u,d t u)\- t= o = (vo,vi). 

Thus, by dJ): 

(19) IMIhi + IMIL" <c(\\cuil{d t uAB)\\ 

Furthermore, the energy inequality on the non-homogeneous Lame system (lemma l2.4f) yields: 



h-v\\ HH(0iT)xn) <C / || curl || i2 dt 
Jo 

||curl (d t u A B)||L-H(o,T) x n) < C ||curl/i(i)||| 2 + ||curl A B)||H- 1((0) r}xn)| . 

which implies, using (|19|l . lemma |2~3|) and the following equation: 

(3d t h + curl curl /i = /3 rot (<9 t ^ A B), 

(20) |!« ||^ + < C ||curlft(t)||i a d* + |K||i= + \\vi\\ 2 H -i | • 

In order to use corollary 12. 1UI we need to add to the left side of inequality l|2UI) the L 2 -norm of h(T). We 
may do so by taking a larger T. Indeed, consider s <E [0, T] such that: 

||Ms)||l== t mm||/i(t)||l 2 . 

Lemma \'2.'6\ gives an a > such that: 

geHDH 1 =► Hfllli^n) < a||curl 5 ||f, 2(n) 

IIM*)lll»(n) < ^/ Hfc(*)ll2.»<n)<ft 



< — 
-TJ 



\\cwc\h(t)\\l 2(Q) dt 



Inequality (|20|l taken with initial time t = s yields: 
E{s) = \ (||u(s)||L HQ) + IIM«)lll»(n) + «ll^)lli"(n) y 

< C I / ||curlft(t)||» a(n) di + |K*)||2, a(n) + ||M«)II 



The energy E being time-decreasing, this implies the inequality of corollary 12.101 and so the uniform 
decay of solutions of 0J. The proof of b) is complete. 

2.4. Polynomial decay. 
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2.4.1. Abstract framework. We shall use here the notations of paragraph s. 3. Tl Let iV be a positive integer, 
and Qt the quadratic form defined by: 

N N 

Qt(z) := £ (ll^z(O)H 2 - ||d'z(T)|| 2 ) = ]T q T (d l t z), D(Q T ) = D(V N ). 

1=0 1=0 

The function d\z being a solution of 1)15(1 . its energy decays with time, so that Qt is positive. Recall the 
definition: (cf T4] ) 

Definition 2.11. The quadratic form Qt is said to be closable when the closure Xq of D(Qt) in X 
for the norm: 

\\*o\\qt = V WMl 2 + Qt{z) 

is complete for this norm. 

Remark 2.12. This is equivalent to the fact that for all Cauchy sequence in D(Qt) for the norm ||.||q t , 
(z fe ), converging to in X, we have: 

lim Q T {z k ) = 0. 

k — > + oo 

We shall again assume the compactness of the embedding: X — > DiJ*)' .. 
The following classical argument goes back to Russel |15| . 

Lemma 2.13. Under the following assumptions: 

a) there exist T, C > such that: 

(21) Vz eD(V N ), ||z(T)|| 2 <C(||^|| 2 _i + Qt(^)); 

b) system U5\) have no non-zero solution of constant energy on [0, +infty[; 

c) the quadratic form Qt is closable. 
There exists C > such that: 

(22) Vzo £ D(V N ), Vi > 0, ||z(t)|| 2 < ^||z ||^. 

Proof. We shall first use a compactness argument similar to that of proposition 12.91 Let T be a large 
positive real number, such that (|21|l holds. Consider Xq, the subspace of X introduced in definition 
(|2.11|) . Extending Qt to Xq by continuity, we can still write inequality (|2T)l for v e Xq. Consider the 
following closed subspace of Xq: 

J T := {z G X%, Qt(z) = 0}. 

By (I21f) we have, for any zq e J T : 

\\M\Qt < C\\z \\-i. 

Using assumption b) as in the proof of proposition l2.9l we obtain that for T large enough: 

(23) J T = {0}. 

From now on, T will be taken such that 12-il) holds. The same process as in the proof of proposition 12. 91 
yields: 

(24) Vz eXg, \\z{T)\\ 2 <CQ T (z). 
Elsewhere, there would exist a sequence (zg) of elements Xq such that: 

(25) ||z fe (T)|| 2 = l, lim Q T (z k ) = 0. 

k — >+oc 

Up to a subsequence, we may assume that (zq) converges weakly to in Xq. By (|23|) and 125|) . zq = 0. 
By compactness: 

— » 0. 
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In view of (|25[) . this contredicts 1|21|) . 
By the triangle inequality, we have: 



\/q T (x + y) < y/qrix) + VQriy)- 
Noting that (Id + V) is an isomorphism from D (V 1 } to X, it is easy to show: 

N 

Qt{z) < CJ2<1t ((H - d t ) l z) 

1=0 
N 

<cj2(\\m\\f-MT)\\f). 

1=0 

From this and J22J, we deduce: 

N 

z G D{V N ) => \\z(T)f < Cj2(\\<0)\\i - MT)\\f) 

1=0 

N 



\In) 



z eX^ \\z(T)f_ N < C£(\\z(0)\\l N - \\z(T)\ 

1=0 

N 

(26) zo e X => \\z(T)f < C\\z(T)\\%J2 (MVWLn - \\z(T)\\In) 

1=0 

The last inequality follows from the trivial bound: ||/|| 2 < ||/||j\t||/||-jv- 
Set: (3 n := £f =0 HnT)]]^. We have: 

Pn< (N + l)\\z(nT)\\ 2 . 

Thus, using Q26|l with the initial data t — nT instead of t = 0, there exists a constant C > such that: 

p 2 n+1 <C\\z(0)\\%(f3 n -f3 n+1 ). 
The following lemma (standard in this setting, see lemma 2.1] completes the proof of lemma 12 . 131 
Lemma 2.14. Let ((3 n )neN be a sequence of positive real numbers and Mq > such that 

P 2 n+1 <M (f3 n -l3 n+1 ). 

Then: 

Vn>l,A,<^. 

n 

Inequality (|26[) implies, with lemma 12. 141 

\\z(nT)f<^\\z(0)fN, 

from which we deduce, taking into account the decay of the energy of z that (|22|1 holds. □ 
Proof of lemma \2.14\ Let a n := nf3 n . Then: 

/ n \ 
n + 1 \^o( n + 1) / 

In particular: 

a n+ i > 2M 4a„> a„+i 
Assume there exists at least one integer n such that a n +i > 2Mo. Let TV be the smallest of these integers. 
Then: 

Oin > Otn+l > 2M , 

which contredicts the minimality of N when N > 1 or the fact that ao is null when N ~ 0. □ 
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2.4.2. Proof of proposition \2. 7| ). Let V be a solution of Q, with initial data in D(A ). The first step 
of the proof is to approach V by a solution U of the Lame system. 

It is easy to see that (d^v,d^ +1 v)^ t =o is in Hq x L 2 . The operator A e being an isomorphism from 
i?o n H 2 to L 2 , the operator C is an isomorphism from D(C) to -ffp x -£ 2 - As a consequence, we may 
choose (uq,ui) such that: 

(uo.Ui) G D(jC n ), (-C) N (u , Ul ) = {d?v,d? +1 v) n = Q . 

The corresponding solution of the Lame system U = (u, d t u) satisfies: 

f d? +1 u n = a = d? +1 v n = a 

Set w = v — u. We will first show: 

N N 

I, 2 ((0,T)xO)- 

1=0 3=0 

We have: 

(28) d 2 w-A e w = K {cm-\h)f\B, (d?w, d? +1 w) = (0, 0). 
Equation implies, for Z G {0, .., TV - 1}: 

(29) ||A e 9^ rt=0 || i2 < \\d l t +2 w lt=0 \\ L 2 + ||curia^ rt=0 || i2 . 
But any g e H 2 (Q) n Ho(fi) satisfies: 

IMIjfi < |M|if= <C||A efl || £ 2, 

which yields, with 129|) : 

(30) Hw n=0 \\ H i < C {\\d\ +2 w^ \\ L2 + ||curia^ rt=0 || i2 ) . 
Since {d^w, d^ +l w) \t=o is null, we deduce from (|3Tj|l : 

w-i 

V/ = 0,...,iV, \\d\ +l w lt =4l, m + \\d l t w tt = \\ 2 H i {Q) <C^ ||^curl/i rt=0 || i2( n) 

N 

(31) V=o|li^) + ||^rt=o||^ J(n) < ||^curl/i|| iW)xn) . 

j=o 

(the second ligne is a consequence of the standard trace theorem with respect to the time variable). With 
the energy estimates of lemma l|2.4|) . applied to l|28[) . we get, for any < I < N: 

II^IIki((o,t)xQ) < c (\\d l t +1 w\t=o\\ 2 L2 {n) + \\d l tW [t =o\\ 2 H t (n) + ||curl^/i||| 2((or)xfi) ) , 

which yields exactly l|27l) . 

On the other side, assumption l|14|) implies 

INI^o) + ||«i||i a(n) < C (||uo||i» ( n) + H«i|lSr-i(n) + Cr(«)) ■ 

Hence: 

\\v f Hi + IMIi> <c (\\uo\\ 2 H i + \\Mh + IKII*» + IMIi») 

(32) <C (\\voWh + INI 2 ^ + Qt («) + M + IKH^ + Kill*) 
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With we get: 

N 

Qt( w ) <c£l|d' +1 Hii 2((0iT)xa) 

N 

<cJ2\\ d3 t cmlh \\hmT)xci)- 

3=0 

Taking into account (|31fl (with I = 0) and the equation: curl (ptvhB) = dth+ -| curl curl h, which implies 
(with lemma l2~3'|l : 



JV 

II »3 1 u\\ 

i 2 ((o,T)xn)> 



Q$(v)<J2\\d!cuTlhr 



3=0 

we deduce from <|3"2"1) : 

IMI»J(n) + ll«illi-(n) < C \ \\vo\\ 2 L 2 (n) + IHIIr-^n) + ll^ curl/l l' 2 



AT 

' ' l '"l 2 ((0,T)x!]) 



3=0 



Now from lemma|H \\h Q \\ 2 L2 < C (||curl/i||| 3((0ir)xn) + \\d t cuilh\\l 2({0jT)xa) 
This gives the following inequality on solutions of Q with initial data in D(A) 




rn^^c \\\v \\i (Ay +^(\\div(o)\\ 

[ 3=0 

It is easy to check, with the criterum given by remark \2. 121 that the quadratic form: 

N 

Qt{V) = Y, / |curl^/i| 2 dy, 



1=0 

with domain D(A N ), is closable. All the assumptions of lemma 12. 131 with V — A hold which completes 
the proof of proposition 12. 71 

3. Defect measures 
Let iV be an integer. For an open subset U of an euclidian space, we set: 

L 2 {U) := L 2 (U, C w ), H S {U) := H S {U, <C N ). 
We consider an open subset fl of R™, n > 1, and a sequence (u k ) of functions on R t x Q, y such that: 

(33) u k — ^ in Hf oc (R x Q), 

k — >+oo 

(in the sense that (ipu k ) converges weakly to in i? 1 (M x fi) for all (p € C^°(K x fi)). We assume that 
every u k is solution of a wave equation in f2: 

(34) (y 2 d 2 t ~ A)u k = 0, in 1x0. 

We shall introduce in this section a measure describing, from a micro-local point of view, the defect 
of compactness in H 1 of the sequence (u k ). This description is of fundamental importance to show 
the observability inequalities of the preceding section, for the Lame system may be decomposed in two 
waves equation (see paragraph 13 . 5 . l|l . Micro-local defect measures have been independently introduced 
by P. Gerard and L. Tatar [HIEj- We shall follow the construction of N. Burq and G. Lebeau, which 
describes the defect of convergence up to the boundary of f2. 

We assume, for the sake of simplicity that the functions u k are smooth, so that their traces on the 
boundary are always defined. In the sequel we shall always reduce to this case. 
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In subsection 13.11 we will give a few definitions and notations. In subsection 13.21 we will state an 
existence theorem of micro- local defect measures and set out their first properties. Subsection 13.31 is 
devoted to the propagation theorem of the measure (proved in 0), and subsection l3.4l to some important 
properties of the traces of u k on the boundary. Finally, in section l3~5l we shall apply the construction of 
the measure to the case of a sequence of solutions of the Lame system. 

3.1. Notations. 

3.1.1. Local coordinates. Consider an open cover of f2: f2 = Ui=o ^j's wnere C Q and, for all j > 1, 
Qj is a small neighbourhood of a point of <9f2, such that on Qj, there are geodesic normal coordinates: 

zen 3 ^(y',x n )eY :=Y'x]0,l[, 

where x n is the distance to the boundary, and Y' an open subset of M™ -1 . Most objects introduced here 
are global objectsm but we will mainly use local coordinates. For a large part of this section we choose 
one of the open set flj, j > 1. 

Set X :=Ix Y and denote the elements of X by: 

X (xq 7 X\ : .., X n — \ , X n ) , Xq — t 1 y — (xi,X2,.., X n — \ , x n ) , 



Let: 



1 "- J — {(x',x n ) e R"+\ x n > 0} , := R^ +1 , dX = X' x {0}, X := X 1 x [0,l[. 



The set X is an open subset of K, . Let g be the natural metric on Y, induced by the change of 
coordinates. In a geodesic system of coordinates, g is of the form: 

" g'(y) o " 
o 1 



g = detg. 



3.1.2. Bundles on X. Let's consider T*X = X x E n+1 the cotangent bundle of X and S*X the spherical 
cotangent bundle, which is defined to be the quotient 



S*X:= (T*X\{\£\=0} J/ ^ + , 

by the action of : (A, ^) i— > A£. The elements of those two bundles will be denoted by: 

p = (x,0, x £ X, ^=K',gel"xl, e = (r,»?). 

There is a natural euclidian norm for the ^-component of T*X: \\ri\\ 2 :— h]g rj. 

We will also consider T*dX := dX x > x R^, the boundary cotangent bundle and S*dX the associated 
spherical bundle. 

3.1.3. Operators in the interior ofQ. Le S™ the set of matrix symboles of degree m with compact support 
in X, which are the functions: 

a(x,£) £ C°°(X x R n+1 , Mjv(C)), 
whose ^-projection is of compact support in X, satisfying the following estimates: 

(35 m ) Va, V/3, 3C aP , \d^a(x,o\ < C af) (1 + \i\T~ m , 

and which have a principal symbol a m (x, ^), homogeneous function of degree m in £, such that a — a m 
satisfies (35 m _i) for large The operator of symbol a, A — a(x, D), is defined by: 

Av{x) := J a(x,0v(0e ix4 <%- 

In order to act on functions which are only defined in X, it is convenient to consider only the set A™ 
consisting of operators A which are of compact support in X , in the sense that A = ipAip for a function 
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if G Cq°(X). An operator in A™ maps a distribution in X to a compactly supported distribution in X. 
We shall denote by cr m (A) the principal symbol of an operator A of degree m. 

3.1.4. Operators near the boundary. Let S 1 ™ be the set of matricial tangential symbols of degree m with 
compact support in X, defined as the functions: 

a(x,£) e C°°(X x R n , Mjv(C)), 

whose ^-projection has compact support in X, satisfying the estimations: 

(36 m ) Va, V/3, 3C Q/3 , |c£cf,a(z, £')| < C«/3 (1 + l£'l) m ~ l/3 ' , 

and which have a principal symbol a m (x, £'), homogeneous of degree m in £' and such that a — a m satisfies 
the inegalities (36 m _i) for large |£'|. We define the operator of symbol a, A = a(x,D'), by: 

M*) : = T^TTT / a(x, OH?, x n )e ix ' 



(27T) 

Here, the Fourier transform of v is only taken with respect to the tangential variable x' . As in the 
interior case, we introduce the set A™ of tangential operators A with compact support in X, i.e such 
that A = tpAip for a compactly supported function tp <G Cq°(X). 

The set of all pseudo-differential operators of interest for us will be denoted by: 

A m := {a = Ai + A b , A, e .4™ A b e -4™} • 

3.1.5. Sobolev spaces. Let s£l and cj be an open set of R™. As mentionned before, we denote by H s (uj) 
the Sobolev space of vector-valued distributions (which may be defined as the set of restrictions to u) of 
elements of H s (W l ), endowed with the quotient norm). We also consider the space H[ oc (u>), the space 
of vector-valued distributions such that: 

V<^ e Cg°(w), ^ueH s , 

and ff^ (w), the space of distributions in H s (oj) compactly supported in ui. The notation H* oc (Z), 
will also be used when Z is not open (Z = R x ft, or Z = X, in the following natural sense: 

u k — ► u ( or = 0(1)) in H? oc (Z) <^=> Vy> e C* °°(Z), <^u fe — ► (= O(l)) in H s , 

where C^°(Z) is the space of C°°, compactly supported functions in Z. We will also consider the following 
spaces, suitable for boundary-value problems: 

H^ c (X) = L\0,hHUX')) 

fflpW = {u e <' C S (A), 3^ e C °°(A), u = ^} . 

Note that the elements of A™ are continuous maps: 

and those of .4™ are continuous maps: 

i_rO,s /^y r \ rrO,s — m (~y\ 

M loc { X ) * "comp ( A j • 

It is possible to "micro-localize" convergence properties in H s and H°' s : 

Definition 3.1. Let p e S*X. The sequence (v k ) is said to be bounded (respectively converging 

to 0) in Hp when there exists A e Af , whose principal symbol is invertible near p and such that (^4u fc ) 
is bounded in H s (respectively converges to in H s . 

Let p' e S*dX. The sequence (v k ) is said to be bounded (respectively converging to 0) in H° p s 
when there exists A £ Af, whose principal symbol is invertible near p' and such that (Av k ) is bounded 
in H°- s (respectively converges to in H°' s ). 
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Note that, according to proposition 17.11 of the appendix, for a sequence of solutions of l|34[l . the 
convergence in H ' 1 and H 1 are equivalent. The spaces H ' 1 and the tangential operators are thus well 
fitted for the description of the H 1 convergence of (u k ). 

3.1.6. Melrose's compressed cotangent bundle. We shall now introduce a bundle which naturally contains 
as subbundles both bundles T*X and T*dX. For this purpose, set T£X := X x R n+1 , endowed with its 
canonical topology and consider: 

t*x t;x 

The mapping j restricts to a continuous map: 

T*X — > T b *X n {z„ > 0}, 

which identifies T*X to a subbundle of dimension 2(n + 1) of the interior of T£X. Furthermore, the 
restriction of j to x n — defines a map from T*X n {a;„ = 0} to T b *X n {x„ = 0}, whose kernel is the 
set {£' = 0}. This clearly identifies: 

T*dX sa (T*X n {x„ = 0})/R c „ , 

(quotient taken by identifiying all the points (i', £',£„), € R) with a 2n-dimensional subbundle of 
T£X. The set of all sections of T£X, with the above identifications, may be seen as the dual bundle of 
the bundle of all vector fields on X tangent to dX. It is called the compressed cotangent bundle. 

We will also consider S£X the spherical bundle of T£X, which naturally contains the spherical bundles 
S*X and S*dX. 

3.1.7. Symbol of P and related manifolds. The equation (|34() takes the following form in local coordinates: 
(37) Pu k = 0, P := -g- 1/2 d Xn g^ 2 d Xn + Q, 

where Q is a scalar tangential differential operator of degree 2. Let q(x, £') be the principal symbol of Q, 
and p(x, £) = + q(x, £') the principal symbol of P. They are both scalar, homogeneous polynomials of 
degree 2 with respect to £. Let: 

Char P := { (X, f ) S T*X, p(x, £) = 0} , Z := j (Char P) , £ := j (Char P) U j ({:£„ = 0}) , 

and S'charP, SZ, SZ the corresponding spherical bundles. 

Decompose T*(dX) (and S*(dX)) into the disjoint union of the elliptic region £, the glancing region 
Q and the hyperbolic region Ti: 

£ := {go > 0}, Q := {qo = 0}, H := {q < 0}. 

3.1.8. Global measure. The defect measure is at first constructed in each of the preceding local coordinate 
systems. The objects obtained are then pieced together to M = lx O. It is easy to define from local 
objects global Sobolev spaces and bundles on M, such as Melrose's compressed cotangent bundle T^M . 
We shall use the same notations (charP, Z, Z, S'charP, SZ, SZ,. . . ) for the local and global objects. 
The definition of global operators is less natural in our setting. The symbol A m will denote the set of 
operators A acting on functions on M , which are of the form: 

,/ 

A = Y, A u)- 

where ^4(o) is a classical pseudo-differential operator of order m with compact support in M and each 

is an operator of the sets A m defined in each system of local coordinates. The global space A m depends 

of the coordinate patches chosen, which shall not cause any problem in the remaining of the article. For 
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a totally intrisic construction, we could have used Melrose's totally characteristic operators (see chap 
18.3]). 

3.2. Existence of the measure. 

3.2.1. The existence theorem. The next elementary proposition shows that for any A G A , the behaviour 
of Au k in H 1 only depends upon the restriction of its principal symbol to SZ: 

Proposition 3.2. Let A b G AT £ . Then: 

A b u k — ► m H l . 

fc^O 

Let Ai G A®, whose principal symbol vanishes on Char P. Then: 

AiU k — > in H 1 . 

k^O 

According to proposition 13.21 it is sufficient to describe the H 1 convergence of (u k ) near SZ, in the 
sense given by definition l3.il Let M. be the set of matice- valued measures on SZ, i.e. the dual space of: 

C : = C ° (SZ,M N (C)) , 

and M. + the subset of all positive measures in A4, i.e. measures (x which satisfy: 

Vz G SZ, b{z) > =^ (/i, b) > 0. 

(M > means M is positive hermitian). 

Before coming to the main theorem of this paragraph, we shall introduce a technical condition on u k : 

Definition 3.3. Let the sequence (u k ) satisfies Ij33f) and Q34[l. We shall say that (u k ) is regular on the 

boundary when one the following equivalent assumptions is satisfied: 

(38a) u k Xn=0 = o(l) in H^{dX), k — > +co 

(386) d Xn u k Xn=0 = o(l) in H^J^dX), k — » +oo. 

Note that this is a very weak condition: the standard trace theorems imply conditions (38) with 0(1) 
instead of o(l). All the sequences (u k ) in this work shall satisfy this condition. For the proof of the 
equivalence between (38a) and (38b) see lemma 2.6]. 

Theorem 3.4. Letu k be such that Iff/jj) . \34\) and (38) hold. Then there exists a subsequence of(u k ), still 
denoted by (u k ), and a measure n G M. + , called micro-local defect measure, such that /i(£ U H = 0) 
and: 

(39) VA# G A\ j G {1,2}, lim (A 2 u k + A 1 D Xn u k , u k ) = (/i, a2+ ^ ai ). 

k— >-{-oo v y f 

In (|39() . the notation (., .) stands for the i 2 -scalar product on M (in local coordinates, it is the scalar 
product on X using the metric g 1 l 2 dydt) and [i is considered as a measure on the subset ScharP of S*X, 
using the canonical map j, which is an homcomorphism: 

S*CharP SZ\H. 

This is made possible by the fact ^l(£ U H) = 0. In the case where (u k ) is not regular on the boundary, 
it is still possible to define \i, but /x(£) is non-null, which makes the statement of condition l|39|l more 
intricate. 
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Remark 3.5. The measure p — pTL{ Xri>0 y may be seen as the standard micro-local defect measure (cf 
[H]) of the bounded sequence (u k ) of Hy oc (M). This interior measure describe the compactness defect of 
(u k ) in Hl oc {M) (in particular, it is null when (u k ) converges to in this space), but not in H^ oc (M): 
p vanishes when (u k ) concentrates on dM, even if it does not converge to in H^ oc (M). On the other 
hand: 

ip € Cq° (M) Lp\W y u k \ 2 dx+ ip\dtu k \ 2 dx — > (p, y) . 



k- 

Thus, fx sees all the Jfj* (M) density of (u k ) at infinity. More precisely, it gives a micro-local description 
of this density: 

p G supp/i <^=> u k — > +0 in H^. 



When p is a boundary point, one should replace H by H° p 

Theorem 13.41 is a new formulation, using lemma 2.7 of of proposition 2.5 of this article. 



3.2.2. A sufficient condition of nullity for p. Let p € SZ an interior point and A G A 2 , whose principal 
symbol is invertible at p. By elementary symbolic calculus on classical operators, it is easy to show, with 
formula 

(40) Au k — > in H~ l => p i supp p. 
The same statement holds in Q: 

Proposition 3.6. Consider an operator of the form: 

A = A D 2 Xn + AiD Xn + A 2 , Aj G A 3 , a,j := cr(A,-), 

such that: 

(41) Au k — ► in H%~ X . 

Assume that (u k ) is regular on the boundary and let p G Q such that a 2 (p) is invertible. Then: 

p (£ supp p. 

Remark 3.7. Proposition 13. 61 is trivial when A — A 2 G A 2 (it is essentially the definition of if-' -1 ). 

Remark 3.8. Near an hyperbolic point, it is more difficult to state proposition 13.61 Indeed it is much 
more relevant to study p, in the set {x n > £o > 0}, near rays in and out of p. (see paragraph U^OJ. 

Remark 3.9. Note that according to the appendix, the convergence to of (Au k ) in the space H^ 1 (X) 
near i would imply (|41l) . Furthermore, the proof of the lemma will show that assumption 141f) is equivalent 
to: 

{JAu k , u k ) — ► 0, VJ G A , with support close enough to p. 

k — >+oo 

Proof. Let j — cr(J). According to l]4Up. pTLt Xn>0 \ is null, near p. The same property remains to be 
proved on p,1r Xn —oy. Let: 

JeA°, ip G C* °°(R) tel que ip(Q) = 1, J e := rp ) J. 
In view of i|41[) and formula 139(1 : 



e 

Letting e goes to 0, the dominated convergence theorem and the fact that is null on the support of 
pl{ Xn = } give: 

(42) (p,TL g J ^-) = 0. 
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Let tp € S® be scalar, positive, and compactly supported near p such that ai is invertible on the support 
of V'j and choose J such that: 

The equality l|42l) then shows that (/x, Ig?/') = 0, which completes the proof using the positivity of p. □ 
3.3. The propagation theorem. 

3.3.1. The generalized bicharacteristic flow. The characteristic curves of the hamiltonian flow of p: 

Hp = d^pd x - d x pd£ 

define a local flow on T*X. The symbol p is homogeneous of degree 2 in £, so that the flow of H p does 
not yield a flow on the quotient space S£X. To get such a flow, we shall replace p by p/r which is 
homogeneous of degee 1. Note that on the support of p, (where r does not vanish), p is null, so that \H V 
and H p / T are equal. Furtermore, the integral curves of ^H p and H p are the same. 

Let S be a small conic open subset of Z = j'(charP). Set go := q\x n =a, gi '■= dx n q\x n =o and 

£° := E n {x n > 0} 

E 1 := n = s n {x„ = 0, g < 0} 

S 2 :=EH {£„ = 0, g - gi ^ 0} 

E fe+3 := £ n {x n = 0, go = gi = ... = <gi = 0, < +1 gi ^ 0}. 

Assume that in S, there is no contact of infinite order between the bicaracteristic curves of P and the 
boundary, which means that for a certain finite integer J: 

(43) 3 J e N, S = (J E J . 

Decompose E 2 in the disjoint union: 

s 2 = g 2 <+ u g 2 '-, g 2 '+ ■.= y 2 n {gi < o}, g 2 '- := s 2 n { gx > o}. 

The set g 2,+ is the set of strictly diffractive points and g 2, ~ the set of strictly gliding points. 
Definition 3.10. Let 7 be a map form a real interval / to S and: 

T(s) = J _1 (7(s)) e 5CharP 

which is defined as long as j(s) ^ H. Such a map j(s) — (x(s),£(s)) is called a ray, or a general 
bicharacteristic curve when 7 is continuous from / to E and for all sq in /: 

• if x n (so) > 0, r is differentiable in sq and: 

r'(ao) = -H p r(s ); 

T 

• if 7(« ) eHug 2 -+, 

3e > 0, Vs e]so - e, so[U]so, so + s[, x n (s) > 0; 

• if 7(so) € 5\^/ 2 ' + , r is well defined and differentiable near so and: 

r'(s ) = -H qo T( So ). 

T 

(Thus, if 7 stays in this region, its spatial projection is a geodesic of the boundary.) 

Under the assumption H43JI . R- Melrose and J. Sjostrand have shown that for any p £ E, there exists 
an unique maximal ray 7 taking values in E such that 7(0) = p (cf chap 24.3]). In the sequel, we 
shall denote by 4>{s,p) the resulting flow (satisfying (f)(0,p) — p). The function p/r being homogeneous 
of degree 1 in £, the flow <fi passes to the quotient and defines a flow on £/R?_. 
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3.3.2. The uniform Lopatinsky conditions. 

Notations. Let S™ be the set of symbols a(x', £') of pseudo-differential operators on dX, with compact 
support in x 1 ', with principal symbol homogenous of degree m in and A™ the set of corresponding com- 
pactly supported pseudo-differential operators (cf paragraph 13. 1 .31 for precise definitions). The Sobolev 
spaces on dX, defined as those on X, shall be denoted by iJJ, Hf oc d , H~ & . 

An approximate pseudo-differential equation on the traces of u k is said to satisfy Lopatinsky conditions 
when it is independent of the equation Pu k = 0. Precisely: 

Definition 3.11. Under the assumptions (|33() and (|34|l . the sequence (u k ) is said to satisfy uniform 
Lopatinsky boundary conditions near p £ S*dX when: 



• if p £ G, 3-B_i € A g 1 such that: 



u 



(44) 



o = £_! (D Xn u k Xn=0 ) + h k 
h k — > 0inffi a ; 

k— y+oo 



• if peH, 3B € A% such that 



61) 



d 

' D Xn u k Xn=Q - Au k Xn=0 = B Q (D Xn u k Xn=0 + Au k Xn=0 ) + h k 
o-(Bq) invertible near p 
h k — > 0inL? a 



A 6 A\, cr(A) = y/q {x',e) = V^T 2 - \\r]'\\ 2 near p. 

( II 7 ?' II 2 = W9'~ lr l' is the natural euclidian norm in the local coordinate system). 

Examples. • The Dirichlet boundary condition, u k x _ Q = 0, or more generally a pseudo-differen- 

tial boundary condition of the form: 



(45) 



u k Xn=0 = B- 1 D Xn u k Xn=Q + h k , B-i S A~ 
h k —* in Hl B , 



where the eigenvalues of cr(B-i) are all pure imaginary numbers near p, is an uniform Lopatinsky 
boundary condition, whether p is glancing or hyperbolic. In the glancing case (|45|l corresponds 
exactly to the definition (|44|l and in the hyperbolic case, both operators Id — A£?_i and Id +AB_i 
are elliptic in p, and it is easy to show l|44l|) . taking Bo — (Id — A_B_i)(Id + Ai?_i) _1 (where 
(...) _1 stands for a parametrix near p). 

• Neumann condition: 

D Xn u k Xn=0 —+ in L% g 

k— *+oo 

is an uniform Lopatinsky condition near any p £ Q. 

• Boundary conditions which are not Lopatinsky conditions in the hyperbolic region are described 
in paragraph 13.4.11 

In the glancing case, a boundary condition of the form 14411 locally implies better estimates than the 
standard ones on the traces of u k , which we shall state in proposition 13.191 This shows in particular 
that {u k ) is regular near the boundary, and that p,(£ ) =0. As a consequence, the bicharacteristic flow is 
defined /^-almost everywhere. 

The set Q 2 ' + is transverse to the bicharacteristic flow. The next result is necessary to the propagation 
of p by the flow, which is treated in the next paragraph. 

Lemma 3.12. Suppose that on S, u k satisfies uniform Lopatinsky boundary conditions. Then: 

p {g 2 -+ n s) = o. 
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3.3.3. The propagation theorem. When an uniform Lopatinsky condition holds, p propagates along the 
integral curves of the bicharacteristic flow. In the hyperbolic region, there is a jump (which depends 
upon the boundary condition). We shall only state a propagation theorem for the support of p, without 
giving a complete description of the propagation of p. 

Theorem 3.13. Let p £ TtDG such that (u k ) satisfies Lopatinsky boundary conditions. Consider a small 
conic open neighbourhood £ of p in SZ such that on S, \44\) ( or ) holds. Then the support of p is, 
in £ invariant by the bicharacteristic flow. 

(cf [2 chap. 3.3, th.l]) 

In other terms, if p € £ is on the support of p, so is the entire bicharacteristic passing through p in E. 
Remark 3.14. Inside M, theorem l3.13l is an easy consequence of the transport equation on p: 

(46) (fi,{p/T,a}), aeC™(Zn{x n >Q}), 

which may be immediately derived, using symbolic calculus, from the elementary property: 

(47) lim (A 1 Pu k - PA lU k , u k ) =0, A x € A\, 

k^+oo 

obtained by integration by parts with the equation 1371) . Near a boundary point, property i|47J) . with 
A\ 6 A\ , still holds with an additional boundary term. Consequently, l|46(l holds only for a certain class 
of function a 6 C^°(SZ), satisfying a particular boundary condition on {x n — 0} (condition chosen to kill 
the boundary terms when k tends to oo). The proof of the propagation theorem, which is fairly technical, 
uses (|46|) . and near strictly diffractive points, lemma 15. 121 The boundary condition on a gives the exact 
value of the jump in the hyperbolic region. See [2j par. 3] for details. 

3.4. Estimates on traces. We now state precise properties of the traces of u k in the hyperbolic, elliptic 
and glancing regions, which are one of the main tools of the proofs of the following sections. Those 
results are fairly classical, and we only shall give a proof (in the appendix) for the glancing case. See 
appendix] for proofs in the hyperbolic and elliptic cases. In this paragraph, we shall always assume (u k ) 
satisfies and 

3.4.1. Hyperbolic region. Near an hyperbolic point, one gains without any boundary condition, half a 
derivative in comparision with the standard traces theorem. 

Proposition 3.15. Let p ETC. Then: 

«r*„=o . = 0(1) m Hlg 

dx n u k Xn=Q = 0(1) m Lj d . 

In view of the propagation theorem in the interior of M, the support of p is, near p, the union 
of incoming rays (integral curves of H p / T along which £ n < 0) and outgoing rays (integral curves of 
H p / T along which > 0). When the sequence satisfies uniform Lopatinsky conditions, theorem 13. 131 is 
equivalent to the fact that if an incoming (respectively outgoing) ray is in the support of p, so is the 
outgoing (respectively incoming) ray passing through the same hyperbolic point. In the opposite case 
where the support of p contains, locally, only incoming (or only outgoing) rays, one gets a boundary 
condition which is in a certain sense "orthogonal" to Lopatinsky uniform conditions: 

Proposition 3.16. Assume that near p € Ti, on the support of p, £„ > 0. Then: 

D Xn u k Xn=Q + Au k Xn=0 = o(l) in Lj d 

AeAl C r 1 (A) = V^HMP 
On the other hand, if near p G TL, on the support of p, < 0, then: 

D x n u k Xn = - Au k Xn=0 = o(l) in Lj 9 . 
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In particular, if \x is null near p, 

D Xn u k Xn=0 — -> in Lj a , u k ! — -> in Hj d . 

3.4.2. In £/ie elliptic region. In £, the equation l|37|l implies a pseudo-differential traces equation on u k : 
Proposition 3.17. Let p G £ and (u k ) satisfy the assumptions of theorem \3.J\ Let M > 0. Then: 

(48) D Xn u k Xri=0 + Zu k Xn=0 QinHf 

fc— *+oo 

E £ A 1 , c i(S) = iy^go = iVll^'ll 2 v2l ~ 2 near P- 

In particular, if a boundary condition independent of l|48l) holds on u k near po (such a condtion is 
called as in the glancing and hyperbolic cases an uniform Lopatinsky condition) , the traces of u k converge 
to in appropriate Sobolev spaces H- 1 . Proposition I3.17| still holds in a much more general case, for 
example if P is replaced by a non-scalar operator P. The principal symbol of S depends again upon 
the principal symbol of P. In the next proposition, we only state a consequence of this fact when (u k ) 
satisfies Dirichlct boundary conditions (which are of uniform Lopatinsky type) . 

Proposition 3.18. Let: 

P := D 2 Xn +Q 1 D Xn +Q 2 , 
where each Qj is a matricial pseudo-differential operator of degree j, with principal symbols q^. Let: 

p(af,£) := q 2 +CnQi 

be the principal symbol ofP, and p = (x',£') be a point of S*dX such that the matrix p(i', 0, £„) 
is invertible for any real number £ n . Consider a sequence (u k ), weakly converging to in H^ oc (X) and 
satisfying: 

Pu k =0, u k Xn=0 = 0. 

Then for all M : 

We shall later apply the preceding proposition in the elliptic zone {£t^£l with notations of paragraph 
I3.5.2f> of the Lame operator <9 t 2 — A e . 

3.4.3. In the glancing region. The strong results of the two preceding paragraphs do not hold in the 
neighbourhood of a glancing point. In this case, one need boundary conditions to get further estimation 
than the standard traces theorem with loss of one half-derivative. In the case of Lopatinsky boundary 
conditions, the results are similar to those of the hyperbolic region. 

Proposition 3.19. a) Let p G Q. Assume that (u k ) satisfies Lopatinsky uniform boundary conditions 
near p. Then: 

(49) < n=0 = 0(1) m fl* 8 , d Xn u% n=0 = 0(1) m L% 9 . 
Furthermore, if fj, vanishes near p, then: 

(50) u k Xn=0 = o(l) m Hlg, d Xn u% n=0 = o(l) m L\ g . 
b) Assume: 

(51) wf*„=o = °( 1 ) m H p,0> 9 x„wf Xn=0 = o(l) in Lj 9 , 

and that every p G Q is not diffractive, in the sense that at least one off the two half-bicharacteristic 
passing through p stays in dfl near p. Then fi — near p. 

Remark 3.20. As seen in propositions 13.151 and 13.1751 point a) holds in the hyperbolic case, where no 
boundary condition is required. 

The proof of proposition 13 . 19l is given in the appendix. 
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3.5. The Lame system. This subsection is devoted to the Lame system with Dirichlet boundary con- 
ditions on an open bounded subset £7 of R 3 : 

{d 2 u - A e u = 0, (t, y) £ R x Q 
uian = 
Uft=o = u £ Hi, dtUf t =o = u 1 £ L 2 . 

In paragraph 13. 5. ll (|52|) is decomposed into two wave equations. In paragraph 13 . 5 . 2l we shall introduce 
the defect measures associated to these equations. Next paragraphs are devoted to a few elementary 
properties of these measures. 

3.5.1. Transversal and longitudinal waves. The natural energy: 

E(t) := I ! {\d t u\ 2 + t i\Vu\ 2 + (\ + fi)\divu\ 2 ) dy 
z Ja 

is time-invariant. Let Eq be its constant value. The nest classical proposition is proved, for example, in 

m 

Proposition 3.21 (Decomposition of the Lame system). There exists a constant C > such that for 
every solution u of \bty . there exists: 

u T £ Hl oc (M) , u L £ H? oc (IT) , 

such that: 

(1) u = Ut + Ul, div ut = 0, curl Mi = 0. 

(2) (df — c^A)ut — 0, where c\ := fi. 

(3) {d 2 ~ c 2 L A)u L = 0, where c\ := A + 2/x. 

(4) For every bounded interval I ofR, of length \I\: 

ll^ll«i(ixn) + l|wr|lHi(/xn) < C\I\E . 

(5) If u is in the space vector generated by a finite number of eigenfunction of C, then: 

MTeC°°(M), u L £C oc (M). 

Definition 3.22. The function ut is called transversal wave, and the function mj, longitudinal. 

Remark 3.23. In the sequel we shall often reduce the longitudinal wave to a scalar function, writing 
Ul = V</2, with: 

cp£H 2 oc (M), |M|H»(/xn) < C\I\E , (d 2 -c 2 L A)^ = Q. 

3.5.2. Measures. Let (u k ) be a sequence of solutions of the Lame system with: 

(uo,u?) — - in Hl(fl) x L 2 (fl). 

k — > + oo 

In view of the continuity of the map introduce in proposition 13 . 2 ll 

(u ,ui) i — ► (u L ,u T ) 

H^xL 2 — » (H\lxn)) 2 , 

where / is a bounded interval, both sequences (uy) and (u\) weakly converge to in H^ oc (M. x O). 
Likewise, the sequence (dt(fi k ) weakly converges to in Jfj^ c (R x fi) (it is more convenient to consider a 
derivative of </? in order to work in a H 1 space, for which the defect measures introduced here are well 
fitted). 

Lemma 3.24. The sequences (uy), (u\), {d t f k ) are regular on the boundary. 
(cf H lemme 4.2]) 
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Notations. Let: 

• fiT> U>h and [i be the defect measures respectively associated (up to a subsequence) to (u^), 
and (dtip k ) by theorem 13.41 

• Ht, Hl, Gt, Gl, £t, £l the hyperbolic, glancing and elliptic region of the transversal and lon- 
gitudinal waves. 

All the calculation shall be carried out in one of the J + 1 local coordinate systems choosen in the 
beginning of this section. We shall make a distinction between the spaces of scalar operators: A m , A™ 
(defined in paragraphs 13.1.31 13.1.41 13.1.81 with N = 1), and the spaces of matricial operators A! n , 
(with N = 3). 

The notation x — (y, i) always refers to local coordinates. When a distinction is necessary, we shall 
write coordinates on fl before the change of variables z = [z\, z%, Z3). This global system of coordinates 
has been chosen so that the magnetic field is vertical: B = (B,0,0). 

Remark 3.25. One may identify /it and hl to measures on S£ (R x Qj with values endomorphism of 
C 3 . 

Remark 3.26. Condition c| 7^ c T means that the intersection of Gt and Gl is empty. 

3.5.3. Link between and/i^. Let x be the local diffeomorphism from global spatial coordinates (zi, Z2, Z3) 
to local coordinates: 

z = x(y)> l x'(y)C = v- 

Proposition 3.27. 

Va e C ° (S* b (R x tt),M 3 (C)) , (fx, = (fi L ,a). 

In particular measures [i and [il have the same supports. 

Proof. Let Aj € A* , j = 1,2, A := A- X D Xn + A . 
Set: I k := {Ad t u k L ,d t u k L ). 
On one hand: 

h = -(d t Ad t u k L ,u k L ) — ► (ml, a)- 

— >+oo 

On the other hand: 

h = -(divA\7d t ip k ,d t <p k ) + o(l) — » 

k — >-t-oo 

The boundary terms of this preceding integration by parts converge to according to lemma 13.241 This 
implies proposition 13 . 2 71 when a is of the form £„a -1 + ao, and then by a density argument for any a. □ 

3.5.4. Polarization of fir and Hl- Let 7r bet the orthogonal projection in C 3 on the line generated by £, 
and 7t_l the orthogonal projection on the plane normal to 

(53) TT(V):=\Cr 2 (X-V)C, 7r ± :=Id C s-7r. 

Projectors tt and tt± are defined by formulas l|53|) . on 

Proposition 3.28. The measure [Il is polarized along the direction of propagation, and fix orthogonally 
to this direction: 

pb L = tthltt, [i T = 7r_L^ T 7r_L. 

Proof. The statement on /xz, is an immediate consequence of proposition ^ .271 To show the statement on 
fj,T> take Aq £ Ao- The nullity of divu T implies: 

o = (A,vdiv4,4) — » ( MT ,^). 

Thus: 

< ht, a7r >= 0, < (Xt, g7t_l >=< Mr, a > ■ 
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The symmetry of \it completes the proof. □ 

Remark 3.29. To get more intrisic formulations of the preceding results, i.e. statements where the two 
coordinate system do not mix, one should have considered u\ and u T as section of the tangent space Tf2, 
and defined measures with values endomorphism of TO, (instead of endomorphism of C 3 ). 

3.5.5. A decoupling lemma. The next result, converting an approximate differential equation on u k into 
two equations on u\ and u T , is of crucial importance in the sequel. As before, (., .) stand for the L 2 
scalar product on 1 x O. 

Lemma 3.30 (decoupling lemma). Let A be a pseudo- differential operator of order 2 of the following 
form: 

2 

(54) a= J2 A ^l: j + A i 

3=~M 

Ai e A 2 , Aj e A 3 . 

Then: 

(55) lim (Au T ,u\) = lim (Au k Ll u^) = Q- 

k — > + oo k — >+oo 

Let A be a (1,3) matrix of pseudo-differential operators, with coefficients of the form j54)) , but with scalar 
operators. Then: 

(56) lim (Au T ,d t (f k ) = 0. 

k — > + oo 

Proof. We shall only prove the convergence to of (Au T ,u'l). The proof of rest of the lemma is very 
much the same. We may obviously assume that the oprators Aj have compact support in one of the local 
coordinate system introduced in paragraph 13. 1.81 In view of the equations: 

-g~ 1/2 d Xn g 1/2 d Xn u k T + Q T u k T = 0, -g-^dx^d^ul + Q L u k = 0, 

where Qt and Ql are tangential, it is sufficient to prove the lemma in the cases j — 1,2 and in the 
interior case. 

First case: A £ A 1 . 

In view of: vj, / i>x it is easy to construct two operators: 

suppcr (*T) = suppo- (*L) n Z T = 0, 

where U is a small open subset of f2 such that there exists a function cp £ C°°(U) satisfying ipAip = A. 
Writing: 

A = Atp = A$ L ip + A^t<P, 
we may assume that the principal symbol of A does not intersect Zt (or does not intersect Zjf). For 
such operators, (|55|l holds trivially. For example, in the first case we have: 

Au k T = 0(1) in L 2 . 

Second case: A £ Ag. 

We know that the support of fJ-T^-{ Xrl =o} is included in Qt and the support of fJ-L^{ Xn =o} m Ql- 
As a consequence, we may write A = A(Qt + @i) where Qt and 0£ are tangential operators of degree 
such that: 

suppct (6t) n Gl = suppcr(e L ) n q t = 0. 

We may thus assume that the support of the principal symbol ai of A is disjoint from one of the two 
glancing sets, say Qt- We have: 

(Au^,u k L ) = ( X (x n /e)Au T ,u k L ) + ((1 - X {x n /e))Au T ,u k L ). 
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where x is a compactly supported function in R equal to 1 near the origin. We first fix e and let k tend 
to oo. The second term of the sum tends to in view of the preceding case. As for the first term, we 
have, u\ being bounded in H l : 

I (x(x n /e)Au^, u k L ) | <C\\x(x n /e)Au^\\ L 2 {QihH -i {x , )) \\u k L \\ L 2 i0dM i (x/)) 
<C\\K'_ l x{x n /e)Au T \\ L 2 + o(l), k -> +oo, 

where A'_ 1 6 Ag , with principal symbol equal to near the support of A. We have: 

IIA' ( U\A fe||2 / (X(x n /s)) 2 |a 2 | 2 \ 

v \ {A k k\\sl (x(x n /e)) 2 \a 2 \ 2 \ 
hmsup \ {Au Tl u L ) \ < (ht, j— m >. 

fc^ + oo x T ||? || ' 

When e goes to 0, the right side of this inequality converges (by the dominated convergence theorem) to: 

(MT^^oiMV^iie'ir 2 ), 

which is null, because Qt and the support of a 2 are disjoint. 
Third case: A — AiD Xn , A\ tangential. 

As in the preceding case, we may assume that <7i(^4i) is disjoint with one of the two glancing sets, say 
Qt- Then: 

(Au T ,u k L ) = (x(x n /e)Au T ,u k L ) + ((1 - X {x n /e))Au T , u k L ) 

The second term converges to zero when fc^oo for the same reasons as in the case A £ A 2 . The first 
term may be written: 

( X (x n /e)A 1 D Xn u T ,u k L ) = (D Xn (x(x n /£)A 1 u T ),u k L ) + (R e u k T ,u k L ) , R e e A 

= ((x(xn/e)Aiu T ), D Xn u k L ) + (R s Ut,u'1) + {boundary terms} 

The boundary terms tend to zero when k tends to infinity because (u^) and {u\) are regular on the 
boundary. The proof may be completed as in the preceding case, letting k go to infinity then e go to 
zero. □ 

4. Sufficient condition 

If 7 is a S-resistant ray defined on a real interval ]a, b[, we shall call life-length the positive quantity 
\t(b) — t(a)\. In this section we use tools of the preceding section to prove the following: 

Proposition 4.1. Let T > 0. Assume that every B-resistant ray in £1 is of life-length strictly less than 
T . Then there exists C > such that for evey solution of the Lame system \bty : 

(57) IKH^i + Kills < C (llcurl^u A B)\\ 2 H - 1{(0 T)xn) + \\u \\ 2 L2 + IM^-i) 

Inequality (|57|l is the sufficient condition for uniform decay stated in point a) of proposition 12. 51 
Proposition 14.1 1 thus completes the proof of the sufficient condition of theorem ^ namely that the non- 
existence of arbitrarily large _B-resistant rays on Q implies the uniform decay of the energy for solutions 
of the system of magnetoelasticity. 

To show (I57f) . we shall argue by contradiction, considering the defect measures /j.t,l of subsection l3.5.2l 
associated to a sequence (u k ) which contredicts (|57J) (cf subsection 14. I|) . The bound on curl(9(U fc A B) 
given by the negation of (|57(l implies a strong condition on the supports of these measures (see subsection 
I4.2|l . In subsection 14.31 we make use of this condition, together with propagation arguments near the 
boundary of f2. Subsection l4.4l completes the proof, using the assumption of non-existence of B-resistant 
rays of life-length larger than T. 
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T)xn) + Ikolli* + 



4.1. Introduction of measures. Assume that Q57H does not hold. Then there exists a sequence (u k ) 
of solutions of the Lame system such that: 

(58) 1 = HugH^ + |K||1 2 > k (||curl(^u fe A B)||| r _ 1((0i 

Up to the extraction of a subsequence, one may assume that (uQ,u k ) weakly converges in Hq x L 2 . 
Inequality l|58l) implies that its weak limit is 0. We may thus introduce the defect measures /it, Hl and 
fi of paragraph 13. 5. 21 associated to the sequences and (d t (p k ). To contredict l|5X|l. we need to 

show that these measures are null. Note that (|5^|l implies: 

(59) curl(d t u k AB) — ► in H~\(Q, T) x Q). 

k — > + oo 

Remark 4.2. By a density argument, it suffices to show (|57|) with (uo, U\) generated by a finite number 
of eigenfunctions of C. We may thus assume, that u\ and u\ are C°° . 

Remark 4.3. We will indeed show a more precise statement than proposition 14. II namely that if (u k ) 
is a sequence of solutions of the Lame system converging weakly to in the energy space and satisfying 
B59JI then the set (supp fir U supp hl) D {t G (0, T)} a an union of £?-resistant rays of length T. 

4.2. Condition on the supports. We may see B as a vector field on f2, i.e. a section of TVt. To avoid 
confusions, the magnetic field considered as a vector field shall be refered as B. In a local coordinate 
system, if: 

X : U C O — > I" 

is the change of coordinates, and x' its differential, B is equal to ^B. Notation B shall always refer to 
the vector of R 3 of coordinates (B, 0, 0). As before, [zi, Z2, Z3) refers to the global spatial coordinates on 
f2, before the change of variable. 

Lemma 4.4. Assume if5^) . Then, on the interval (0, T), f/ie transversal measure charges set of all points 
whose direction of propagation is orthogonal to B and the longitudinal measure charges set of all points 
whose direction of propagation is parallel to B. 

(60) /i T ll(o,T) =Mt1 (0 ,t)%x, B 1 - :={(t,y,T,r]); f Br 1 = 0} 

(61) Ml1(0,t) = j«il(o,T)l^//, B /f := {(t,y,T,n); n E vect(gB)}. 
Proof. Set: 

_^ / -d Z2 u 2 - d Z3 u 3 

(62) i?u := curl(<9 t u AB)=9 ( <9 2l u 2 

Transversal measure . The measure /iT does not charge neither TCt nor £t- Thus, it suffices to check: 

supp/iTl(o,T)l{ :E „>o}ueT c B± . 
Near the boundary, by proposition 17. II of the appendix, 

Ru k — > 0inL 2 ([Q,l[,H^ c (X')). 

k — > + oc 



Thus, according to the decoupling lemma, whether Aq has support in the interior or near the boundary: 
(63) VA) e A , suppAo c{i£ (0,T)}, (A i?w T ,u T ) — ► 0. 

Because divwj. = 0, formula Q62)) may be written: = d t d Zl u^. This implies (by paragraph] 



(0,T)}n{(7 2 (9A l ) ^0}) =0. 
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This clearly shows the announced result near an interior point. When p G Qx one may write: d Zl — 
fod x „ + Fit where fo is a function and F\ a first order tangential differential operator, and the following 
basic fact completes the proof: 

a(F 1 )(p) = ^ fj'±B. 
Longitudinal measure. The first coordinate of Ru L is: 

-d t (dl+d 2 Z3 W k . 

Its scalar product with dtf k gives: 

VA e A , supp (A ) c{te (0, T)} => lim (A (d 2 Z2 + d 2 z )d t ^ k , d t ip k ) = 0. 

k — > + oo 

This implies (again by paragraph l3.2.2ll . that \i (thus ^l) vanishes, in (0, T), on the set of all p such that: 

p G {x n > 0} U Q T , a 2 (d 2 2 + d% )(p)?0. 
Hence JSJ. □ 

4.3. Support of the measure near the boundary. For any symbol go with support in {x n > 0}, we 

have: 

(fir, H PT / T q ) = (p L ,H pL / T q ) = 0, 

which shows, in the interior of f2, the invariance of each measure by the appropriate hamiltonian flow. 
Unfortunately, the condition: 

(64) u Tldn +u k Lldn = 

is not always sufficient to extend such a property in the neighbourhood of dfl. Indeed, without any 
additional assumption, /i£ and /it are not deterministic: the value of the two measures for time t < to is 
not uniquely determined by their value for time t < tg. In our case, this convenient property holds thanks 
to the strong conditions on the support of px and fiL- As announced before, we shall only describe the 
propagation of the supports of the measure. 

Lemma 4.5. Let px and ph be the defect measures associated to a sequence of solutions of the Lame 
system satisfying ^5iJ\) . Let: 

p=(x',i') = (t,y',f,fj')eS*dX 
and n the unitary exterior normal vector to dfl at y' . Them px and pl both vanish near p except possibly 
in the following cases (cf figure^): 

(1) pl is null. The support of px propagates along the transversal flow and: 

• (Hrw) case: p £ Tlx, fj' = and B is orthogonal to n; 

• (Ttxi?]) case: p G Ttx, fj' 7^ and B is normal to the reflection plane; 

• (Gxw) case: p is diffractive for the transversal wave (i.e. p G Qx and the bicharacteristic 
ray passing through p only intersect the boundary at p), and B is orthogonal to fj' ; 

• (£7t[2]) case: p G Qx is not diffractive for the transversal wave, and B is normal to the 
reflection plane. 

(2) px is null, the support of px propagates along the longitudinal flow and: 

• {TLl) case p G TLl, fj' = and B is parallel to n; 

• (Gl) case p is a diffractive point for the longitudinal wave and B is parallel to fj' . 

(3) Both measures px and px are non null, p G 7ix H TLl and: 

• (T — > L) case: B is orthogonal to the transversal ray coming in, and parallel to the longitu- 
dinal ray going out of p. The support of px is an union of incoming transversal rays. The 
support of px is the union of all outgoing longitudinal rays going out of points of Tlx ^TLx 
where the transversal rays of the support of px come in; 
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Figure 2. Cases arising in lemma ET51 



• (L — > T) case: B is orthogonal to the transversal ray going out of and parallel to the longi- 
tudinal ray coming in p. The support of fiL is an union of incoming longitudinal rays. The 
support of fix is the union of all outgoing transversal rays going out of points of TLt C)TLl 
where the longitudinal rays of the support of p,L come in. 

All the assertions of lemma ET51 should be understood in a neighbourhood of p. The reflection plane at 
a boundary point p = (t, y', r, if), defined as long as rj' ^ 0, is the plane passing through y' and generated 
by n and rj' , thus containing the bicharacteristic ray passing through p. The statement 11 B is parallel to 
Tf'" must be understood as "the vector gB of the cotangent bundle of f2 is parallel to 77"'. 

Notation. Let p be an hyperbolic point for the transversal (respectively longitudinal) wave. We shall 
denote by £^ (respectively the incoming and outgoing vectors through p: 

& := ( Ut = vk^¥f¥ ) ir = { -L ) 



(respectively with "L" instead of "T" ) . 
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We shall write 77^, rj^ the spatial components of this vectors. For example: 

Let's postpone the proof of lemma FOl to show, as stated in the introduction of this article, that in 
the (T — > L) and [L — ► T) cases, the angles of refraction and incidence have a fixed value, determined 
by the quotient ct/cl- Consider for example the (T — * L) case. Let ar be the angle of incidence of the 
transversal wave, /3l the angle of refraction of the longitudinal wave, or and cll the following numbers: 

ay := tanc^T = ,, , o-l := tanpx 



WW WW 

The incident and refracted waves are orthogonal, so that: 

(65) ||rf || 2 - i nT i nL = 0, i.e. a T a L = 1. 
Furthermore, the definition of £„x and £„t yields 

4\m 2 -r 2 = o, cira 2 -T 2 =o, 

which gives the equation: 

(66) 4(l + 4)=4(l + ai). 
Equations l(65jl and Q66J) imply the formula announced in the introduction: 

olt — arctan — , @l = arctan — . 

ct c L 

By a similar calculation, one gets, in the (T — * L) case: 

0.1, = arctan — , /3t = arctan — . 

CL CT 

There are thus very strong constraints for the possible transfer of energy from one wave equation to the 
other. 

Proof of lemma \J~5[ Case (J3J: fix 7^ 0, fiL 7^ near p. 

In this case, p ^ St U £l- It is also easy to show that p £ Gt ^ Gl- Indeed, if p £ Gt then it also 
belongs to TIl (it cannot be a point of £l, and Gl and Gt are disjoint). But pt being non-null near p, 
fj is orthogonal to B (by lemma 14.4(1 so neither ff^ nor 77^ are parallel to B , which implies (again by 
lemma l4~4"ll that p^ — near p, contradicting our assumptions. Likewise, if p G Gl, fj' must be parallel 
to B and pt null near p. Thus p € 7Yt Pi TIl- The support of measures /jt and pl is, near /5, an union 
of incoming and outgoing maximal rays. 

Let's first assume that the support of pl contains the ray going out of p. Then: 

so that 77^ is not orthogonal to B . As a consequence, the support of pt is only made of incoming rays, 
and the fact that pt ^ implies: 

i] T ±B. 

Thus rfj^ is not parallel to B . This is the (T — > L) case, and it remains to show the statement of the 
lemma about the transfer from transversal incoming waves to longitudinal outgoing waves, which may 
be formulated as follow: for any p £ TLt H TLl near p the following equivalence holds: 

(67) /jesupp/iT /3esupp/j£. 
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Let's assume for example p ^ supp px- Then px is null near rays coming in and going out of p and by 
the hyperbolic theory (see proposition 13 . 16fl : 

u L\x n =Q = - u T\x n =0 , r 4 m L? d . 

By the propagation theorem, the support of ph propagates near p. But this support is an union of outgoing 
rays. Consequently, it is empty near p and px is null near p. The implication p (£ supp px p (f. supper 
may be shown in the same manner. 

If the support of p^ contains no outgoing ray near p, it must contain incoming rays. This corresponds 
to the (L — > T) case, which may be treated as the (T — > L) case. 

The study of all other cases relies on a transfer argument on boundary conditions, stated in the 
following technical lemma [OH rougly, a boundary condition on the longitudinal wave implies one on the 
transversal wave and vice versa. 

Notation. Let v TL be the functions u T L considered as vector fields on fl. In local coordinates, if \ 
denotes the change of coordinates, we have: 

„fe \ / „,fe 



v 



Tl \ / LI 



'Tn J \ V L 

Lemma 4.6. Let p G S*dX, and (u k ) be any sequence of solutions of the Lame system weakly converging 
to in Hl oc (R x H). 

• Assume the following aproximate equation for some A\ G A 1 : 

(68) d Xn v Tn[Xn=0 = A lV * nlXn=0 + o(l) m L% 9 . 
Then: 

EE1) &v'V% n =o = -Aid Xn <p% n=0 + o(l) tn Hl d . 

• Conversely, if, for some A_i G A^ 1 the following equation holds: 

(69) ^ Xn=0 = A^8 Xn ^ Xn=0 + o(l) m H 2 p d . 
Then: 

ESI) Q«„«r»rx„=o = - A V^i"Tn tl „ = o + o(l) m L% 9 . 

• Moreover if in addition to \69ji . fj ^ and a(A-i)(p) ^ 0, then: 



GST) 4rx„=o = Z-i3x„4rx n =o + o(l) /fi 



/ . 
Z_i G AT 1 , er_i(Z_i) = ||r/||~ 2 %9 |_ m 

V a(A^)- 1 

Proof. First note that the Dirichlet condition on u k implies: 

(70) d*„4nr*„=o = &v'V% n =o + 0(1) in H^(dX), 

Indeed, the equation div u T — implies, in local coordinates: 

9x„4„ + + 9 y A 2 = 0(1) in H^(X). 

Which may be written, using tr _ = 0: 

9x n 4„rx„=0 = fyi«*ir*„=o + ^fii»L2r«„=o + 0(1) in H^(dX), 
yielding i|70[l by the definition of ip. 
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Assume (JBSJl. By (J7UJ) and the nullity of v^ Xn=Q : 

Atftf Wn =o = -Aiv k LnlXn=Q + o(l) = -A x d Xn ^ Xn=0 + o(l) in L\ d . 
Now assume Hence: 

Which implies (|69ll) using l|70() on the left side of the equation, and the Dirichlet condition on v!f n on its 
right side. 

If, in addition to the assumption (|69[1 . fj' and a{A_ 1 )(p') are non zero, both operators A_i and A y i 
are elliptic at /5, and equations i|69t[l and l|7U|) may be rewritten: 

(71) 4„r*„=o = ^-i^n4nrx„=o + o(l) in L% 9 a(F_i) = UrflU-V^!)- 1 , 
^ n=0 = ^- 2 9x n 4„rx„=0 + 0(1) in Hl'l a(E_ 2 ) = -\\rf\\- 2 near p 

(72) ( 1 ir -=° ) = 9'- 1 ( ^4«=° ) = Z^t^o + 0(1) in Hl'l = . 

<7(Z_ 1 ) = -i||r / '||-V- 1 ( £ ) near p. 
Assertion (|69l'[l is an easy consequence of <|71[1 and (I72[l . □ 

We may now study cases Q and (J2J) of lemma 14.51 

case Q: assume /ij, = 0, p,T ^ near p. There are three possibilities: 

• If /? £E 7ii,, the nullity of /iz, implies, by standard hyperbolic theory: 

u T\x n =Q = ~ u L\x n =0 , i n Hp,d- 

k— >+oo 

So near p, the support of /iy propagates. It is easy to see that condition l|60() on the support of 
Ht implies, if pr does not vanish, that this is one of the four cases described in the JQl of lemma 

ED 

• Assume p £ £l. The standard elliptic theory ( proposition 13 . 1 7|l implies: 

d Xn ^ Xn=0 + &p* Xn =o — in Hl lXn=0 . 
With lemma FOfl this yields the following equation on the traces of u^: 

(73) 4rx„=o = iZ-iD Xn u^ Xn=0 + o(l) in H^, Z_i = x'^-xx' 

where the principal symbol of the operator Z_i S «4.~ is given by (EMJ. Notice that the 
eigenvalues of tr_i (iZ^i), thus those of rr_i(iZ_i) are pure imaginary numbers. As a consequence, 
the boundary condition (|T3|> is an uniform Lopatinsky boundary condition near p (see the example 
following dcfmition l3.11l) . which shows again the propagation of px- As in the case where p £ Hl, 
it is easy to see that this is one of the four cases of lemma 14.51 

• The case p € Gl is the most difficult. When pr is non-null p, must be in Ht- We use a 
contradiction argument to prove the propagation of the support of pr- Let p 6 H.t such that 
the ray coming in p is in the support of px, but not the ray going out of p. According to the 
standard hyperbolic theory ( proposition 13 . 1 B)l : 



D Xn u T[Xn=0 - A T 4 r —> Oini? a , cti(A t ) = ^v t t 2 - ||r/|| 2 . 
This implies, by lemma FTBl a boundary equation on dt(p \ of the following form: 
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which is an uniform Lopatinsky condition near p because p E Gl- In view of proposition 13 . 1 9l on 
traces in the glancing region, such an equation implies, with the nullity of p near p the following 
conditions: 

dt<P\ Xn =o in H^g, d n d t <p k Xn=0 -> in L\ a . 
The operator dt being elliptic at p, this shows that u k L tends to in g , and thus: 

u T\x n =Q , ~* in Hj a , 

Hence the propagation of the support of px near p, which contredicts the assumption on rays 
coming in and going out of p. 

Similar arguments show that if the ray going out of p is in the support of px, so is the ray 
coming in p. This proves that the support of px propagates near p. Notice that this is necessarily 
the (Htm) case. 

Remark 4.7. Case JJJ, which appears in the study of linear thermoelasticity, was precisely described in 
EJ. The authors show a result of propagation of px, determining all the characteristic elements of this 
propagation, which gives in particular the polarization properties of px- In the case of the system of 
magnetoelasticity, the polarization causes no problem by and it suffices to show the propagation of the 
support of px (or that of px in case (2)). As mentionned in the introduction, one may consider that 
the only component of px and px which is resistant to the dissipation is the component parallel to B, 
cancelling the quantity uAB. 

Case El we assume now that px — and px ^ 0. We argue in a similar way, considering three possibilities: 

• If p E Ht, the standard hyperbolic theory gives an approximate boundary equation on u T , which 
implies: 

"£r*„=o t — > OinJTig. 

As a consequence, the support of px propagates. On this support, ?y is parallel to B, which 
shows, as stated in lemma ET31 that rf = and p E Ttx, or 77'// B and p E Gl- 

• If p € Ex-, we write (as in the similar situation when p £ £■£), the boundary equation of the elliptic 
region: 

D Xn u T[ + E T u T[ — > in Lj a . 

This implies in view of lemma an uniform Lopatinsky boundary equation on dt<p k , thus the 
propagation of the support of p, which is the same as that of px- The fact that 77' 7^ shows that 
p cannot be hyperbolic for the longitudinal wave (in this case outgoing and incoming directions 
are not parallel, thus at least one is not parallel to B). Consequently, p E Gl- More precisely, 
it is a diffractive point: the bicharacteristic passing through p must stay parallel to B , thus its 
direction is constant which is not possible for gliding rays because has no contact of infinite 
order with its tangents. We are in the (Gl) case of lemma 1431 

• If p € Gt> then p € Hl- The fact that p € Gx implies that fj' ^ , so directions r)\ and fj~[ 
cannot be both parallel to B . Consequently, the support of pl is an union of only ingoing rays 
(or only outgoing rays). This gives a boundary equation of the following form: 



d x n d t ^p k Xn ^ + iA L d t tp k Xn=0 — ► in L 2 



p,<9' 



where 1 E {+1,-1}. Notice that dt is elliptic at p, so that we may rewrite this last property 
taking out all the dt and with H 1 instead of L 2 . This yields, in view of lemma EHJ1 a unform 
Lopatinsky boundary condition on u T . The nullity of px gives as before (by proposition 13 . 1 9|l : 



fe— H-00 
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so that propagates, and in view of the particular form of its support, vanishes near p. This 
shows that this particular situation (p £ Qt and fi^ 0) is impossible, and completes the proof 
of lemma 14.51 

□ 

Definition 4.8. We shall call S-admissible points the point of the boundary of S^M which are of one 
of the eight types described in lemma |4~B1 

4.4. Conclusion of the proof. Let S := suppp^ U supper, and BR the subset of S^M, of all points p 
satisfying one of the following properties: 

• x n > 0, p G SZt and r]//B; 

• x n > 0, p 6 SZl and rj-LB; 

• x n — and p is _B-admissible. 

Let $t(p, s) and &l(p, s) be the bicharacteristic flows for the transversal and longitudinal waves. We 
shall define a local continuous flow on BR, denoted by: 

in the following way: 

• if $ In ( l o, s) > and $(/), s) € SZt, or if s) is a B-admissible boundary point of the form 

of lemma |4~51 $ is near (p, s) the restriction to BR of the transversal bicharacteristic flow; 

• if $ Xn (p, s) > and Q(p, s) £ SZt,, or if <&(p, s) is a £?-admissible boundary point of the form 
J2J of lemma |4~51 $ is near (p, s) the restriction to BR of the longitudinal bicharacteristic flow; 

• if s) is .B-admissible of type (T — > L), then: 

$(p, r) = $ T (p,r), if r < s 
$(/9,r) = $ L (p,r), if r > s; 

• if <£>(p, s) is .B-admissible of type (L — > T), then: 

$(p,r) = ® L (p,r), if r < s 
r) = 3>t(p, r), if r > s. 

In view of lemma FOl S (which a subset of BR) is stable under the flow $ on (0, T). Furthcmore, if for 
some p G BR, 

^(P> s ) — I P & BR, SG (0,T), 

s — >s 

(thus p is a boundary point which is not .B-admissible), then p is not in S. Consequently, S is an union 
of B-resistant rays of life-length T. The assumption of non-existence of such rays made in proposition 
14. II shows that S is empty, which completes the proof. 

5. Necessary condition 

Proposition 5.1. Assume that for all T > 0, there exists a B-resistant ray of life-length T . Then for 
all T > 0, there exists a sequence (u k ) of solutions of the Lame system such that: 

(74) l|fit«fUolli» + lKt=ollL — » 1 

(75) \\u k AB\\ H i mT )xQ) , — ► 

Corollary 5.2. Under the assumption of vror)osition YE~i\ the energy of the solutions of the magnetoelas- 
ticity equations does not decay uniformly. In other terms, the necessary condition of theorem^ holds. 

Corollary 15. 21 is a direct consequence of point b) of proposition 12. 51 
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Proof of vrovosition Tb~~l\. This proof is very much inspired by that of theorem 4 of Denote by Z one 
of the indices T or L and set Pz '■= A — v z d^. We start by an elementary remark: 

Remark 5.3. If E is a vector subspace of C 3 and tte the orthogonal projection on E, the defect measure 
of iTEU k z is tt e /i z tte- Furthermore, P z tteu z — 0, so that theorem of propagation 13.131 holds for the 
measure He^z^e if an uniform Lopatinsky boundary condition holds on iteu z . Notice that any scalar 
uniform Lopatinsky condition on u k z yields such a condition on iteu z . If tt e /j, z tte — Hz, the measure /iz 
will be said to be polarized along E. If E is the line generated by a vector H of C 3 we shall also use 
the phrase "polarized along H" . If both measures fj,T^]-e,T+e[ an d ^l1]-£.t+s[ are polarized along B, 
then condition H75|l is fullfilled. 

Let X" > T. Consider a B-resistant ray defined on an open interval / of length T": 

7 (s) = (t 1 (s),y 1 (s),T 1 {s),rj 1 {s)) = (z 7 (s),£ 7 (s)) 
If T" is large enough, then one of the two following assertions holds: 

a) 7(1) contains an interior point; 

b) y~f{I) = r C dfl where F is a closed curved, contained in a plane P which is normal to _B, 
boundary of a convex subset of P, and such that on T, n is orthogonal to B. 

Case b) occures when there exists an infinite boundary B-resistant ray. This case reduces to case a), 
choosing a transversal ray contained in P which only meets the boundary at hyperbolic points. 

Thus, we may assume that has an interior point. We may also assume, possibly moving the origin 
of coordinates, that this interior point is 7(0), and that (i 7 (0), y 7 (0)) = (0,0). Recall that the magnetic 
field is vertical: B = (B, 0, 0). We shall denote by -T~ and T+ the extremal points of/: / = (-T~,T + ). 

If ?y 7 (0) is parallel to B (i.e. if 7(0) is in the longitudinal characteristic set), choose a non-zero function 
ip € Qj°(Q), and set: 

ip k (y) = K^k-s^e^tpiVky) 
ito=V(y5 fc , u k =ikc L UQ. 
Where u k is the solution of the Lame system with initial data: 

{u k ,d t u k ) n=0 = (ug.uj). 

Then: 

llugiiffi — ► K^yw^, \\u k \\ L 2 — ► J ft:- 1 ||^|| i 2. 

k — >+oo k — > + oo 

Thus, condition (|74fl is fullfilled with an appropriate choice of K. 

For small t, by finite speed of propagation for the wave equation, u\ has compact support in Vl. Thus 
«y = and u'l = u k . As a consequence, for small t: 

(1) (J-t = 0; 

(2) fiL is polarized along B; 

(3) the projection of the support of fi^ on M t x £1 is contained in ar 7 (7). 

If ?y 7 (0) is orthogonal to B, we construct a sequence of solution of the Lame system, with the following 
initial data: (cf 2 ) 

Uq = curl (0,0, -ip k ), u k ^ikc T UQ. 

In this case, condition l|74|l is fullfilled for an appropriate K and the defect measures satisfy the following 
properties for small t: 

(1) ML=0; 

(2) [It is polarized along B, 

(3) the projection of the support of fix on R t x is contained in x 1 {I). 
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To show (|75f) . we shall prove that both measures px and pl, are, for t £ J, polarized along B. For < > 0, 
we shall denote by V(t) the following property: in a neighbourhood of [0,i], both measures /i^ and 
/zt are polarized along B and the projections of their support on R t x f2 are contained in 

x 7 (I). 

Let T be the set of t in [0,T + ) such that V(t) holds. By its definition, T is an open subset of [0,T + ). 
We have just shown that in T. We shall now prove that T is closed. We shall use the next lemma: 

Lemma 5.4. Let p = (t,y, f, fj) £ S£M. If px^ t <i and px^ t< i vanish in a neighbourhood of p, so do 
both measures px and px- 

This is a trivial assertion in the interior of fl by the propagation of both measures. Near a point of the 
boundary of tt, one may show lemma I5T41 using the Dirichlet boundary condition on u k and the theorem 
of propagation 13 . 1 31 together with the same type of arguments as in lemma l4~^l 

Let sq > such that V(sq) holds for s < sq. We must check that 'P(so) holds. Three cases arise, 
depending on the nature of p :— 7(so)- 

i) p is an interior point. 

"P(so) is obvious in view of the propagation of both measures in the interior of f2. 

ii) p is of the type Q of lemma 14.51 

This case, where pl vanishes for time t < i 7 (so) near i 7 (so), were studied in [2]. The authors show 
that pl remains null for times greater than £ 7 (so) and that px propagates near 7(so), in such a way that 
in our case, its polarization along B is preserved. In particular property 'P(so) holds. 

iii) p is of the type (Hl) of lemma 1431 y 7 (so) € dfl, ry 7 (so) = and n//B. 

In view of lemma l5~4l the support of the measure pl is contained, near p in the union of the longitudinal 
ray coming in p and the ray going out of p. The support of px, if not empty, if the transversal ray going 
out of p. Let E be the plane orthogonal to B in C 3 . The polarization of px along B shows that 
-k* e pl^e = and thus, by remark f5. 31 and the standard hyperbolic theory of proposition 13.151 

n E UTf Xn =o = - n EU k LlXn=o —>0 in H^ g 
which implies, using again remark I5"3l that tt* e pt^e propagates along the transversal flow near j(sq). 
Thus tt e pt^e vanishes near p. But px is polarized orthogonally to its direction of propagation which 
is exactly B on the support of px near p. This show that px vanishes near p, completing the proof of 

iv) p is of the type (Gl) of lemma 1431 p is a diffractive point for the longitudinal wave, and ?7 7 (so) is 
parallel to B. 

Then p £ Ti.x U Ex- Furthermore, in the longitudinal hyperbolic case, px^t<t 1 (s ) vanishes near p. 
Thus, according to standard elliptic or hyperbolic theory (cf propositions I3.1b1 and I3.l7fl . u T satisfies a 
boundary condition of the following form: 

A.„«rr»„=o = ^«rr»„=o + (!) in L %d, AeAl 

o~i(A) — — i\J\\if\\ 2 — v?pT 2 in the elliptic case, 
o~\(A) = —ij v?pT 2 — \\rj'\\ 2 in the hyperbolic case. 
Each of this equation yields, in view of lemma |4~B1 an uniform Lopatinsky boundary equation on (p k : 
(76) </>f x „ =0 = 5-i A^t„=o + "(I) in H ld- 

As a consequence, the support of p (and that of pl) propagates near p. The polarization of ph along B 
is immediate. The nullity of px near p remains to be checked. This is a general property in the elliptic 
case p £ Ex- In the hyperbolic case, first note that the propagation theorem of Burq and Lebeau 
th. 1] implies with boundary condition (|76J) that p is invariant by the longitudinal flow near diffractive 
points. So the total mass of px is preserved by time, for t close enough to t 7 (so). The next lemma, which 
is a measure version of the conservation of energy for the Lame system, completes the proof of 'P(so): 
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Lemma 5.5. Let if e Cg°(R). T/ien: 

In oi/ier terms, the total mass of the measure (pr + /J>L)\t=s is well defined, and does not depend on s. 
Proof. 

(d 2 - A e )u k = 
Re J d 2 u k u k tp(t)dx - Re J A e u k u k : <p(t)dx = 0. 

Set: V e u := (/xVu, (A + /i)divu) 6 C 4 . A simple integration by parts yields: 

/" y/(t) |^u fc | 2 da: + J ip'{t)\V e u k \ 2 dx. 
Using another integration by parts, and then the decoupling lemma ESDI 

d 2 u k u k ip'(t)dx - / A e u k u k (p'(t)dx = o(l) as k — > +oo 



(<9 t 2 + 4A)u^vj'(i)dx + / (<9 t 2 + c%A)u k L u k L v'{t)dx = 



When fc tends to oo, we get, by the definition of pt and /ij,: 

( w ,^H!».'W> + ( w ,^*/W> = o. 



This completes the proof, noting that on the support of pr (respectively /it), Cr||?j|| (respectively CL||?y||) 
is equal to r. □ 

Lemma l5.5l and the mass conservation of /i£ as time goes by imply that the mass of px is also preserved 
near p, which shows that px vanishes in a neighbourhood of p. 
v) p is of the type (L — > T) of lemma PTBI 

In view of lemma 15^1 and of the assumption P{s) for s < so, the support of pl is contained in the 
two longitudinal half-rays passing through p, and that of pt is only contained in the ray going out of p. 
To prove P(sq), it remains to show that (1l = along the longitudinal ray going out of p. We shall do 
so by a simple polarization argument. Let H (respectively J) be an unitary vector of C 3 parallel to the 
direction of the transversal (respectively longitudinal) ray going out of p. The polarization of pt shows 
that tt^pt^h is null, so that tt^pl^h propagates near p. Now, H is orthogonal to the direction of the 
longitudinal ray coming in p, so that tt^pl^h vanishes, along incoming rays but also, in view of the 
propagation, along outgoing rays. Furthermore pil t>t (s ) i s polarized along J. It is easy to see that 
this last measure vanishes. Indeed, the polarization of p^ implies: 

lt>t T (so) 7r J 7r ff 7r jMi 7r J 7, 'H 7r ./ = lt>t 7 (so) 7r J 7r HMi 7r -H 7r ./ = 

But: 

< H.J > 2 

*J***J = |ff|2|J|2 *J 

Noting that H and J are not orthogonal this yields the nullity of piXt>t (s ) m a neighbourhood of p. 
vi) p is of the type (T — > L) of lemma 

One may argue as before, showing that for every vector K orthogonal to J, tt^ptttk — near p, 
which implies the nullity of pT^t>t^ {s()) near p. 

The proof is completed by reversing time, which yields V(s) for — - T~ < t < 0. □ 
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case v) 



Figure 3. Cases iii), iv) and v) 

6. Polynomial decay 
We shall prove in this section the following proposition: 
Proposition 6.1. There exist T, C > such that for every solution of the Lame system with initial data: 

{u Q ,ui) e D{C N ), 

we have: 

(77) \\uo\\h& + Kill* < c (Ct(«) + IMl£' + H«iHir-x) , 

where: 

• N = 1 if there does not exist on f2 any boundary B -resistant ray of infinite life-length; 

• N = K elsewhere. 

The integer K > 2 was introduced in theorem 

K := sup KfTj), 

3=1. .m 

where ri,...,r m are the spatial images of infinite boundary B-resistant rays and N(Pj) is the minimal 
order of contact of dCl with its tangents parallel to B at points of Tj . The quadratic form Qjj is defined 
in section |21 by: 

Qt( u ) ^J2\\ cml ^ +luAB ^H-H^T),ny 

1=0 

Inequality l|77|l is precisely the sufficient condition of polynomial decay given by proposition 12.71 which 
completes the proof of theorem [21 As in section we shall argue by contradiction, using the defect 
measures of section |31 
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6.1. Introduction of measures. First note that by a density argument, it suffices to show l|77l) for 
initial generated by a finite number of eigenfunctions of £. Assume that (|77|l does not hold. This yields 
a sequence u k of smooth solutions of the Lame system such that: 

(78) 1 = + H\\h > k (Q$(u k ) + \\u k Q \\h + 

As in section 0J we may assume that u k converges weakly to in H^ oc (M. x fi) and introduce the defect 
mesures of subsection 13.51 Note that (|7St implies: 

(79) <#(«*) — > 0. 

k — >+oo 

As a consequence, condition (|59|l of section 0| is fulfilled. In this section, we proved (cf remark rO|) that 
this condition implies that in the interval (0, T), the supports of fix and pl are unions of B-resistant 
rays. 

We would like to show, as in section^ that both measures \xl and /it are null, which would contredict 
(|78J) . We shall first prove by invariance arguments that, as long as N > 1, both measures are null in 
the interior of £1 (subsection 16. 2|) , which yields proposition 16.11 in the favorable case of non-existence of 
boundary B-resistant rays of infinite life length. Subsection l6.3l studies those boundary rays, using traces 
theorem to restrict <|79ll to the boundary of Q. 

In all this section, z — (z\, Z2, 23) denotes a global spatial orthonormal coordinate system of R 3 , in 
which B = (B, 0, 0). As before, notations y and x = (t, y) shall only be used for local coordinates. 

6.2. Nullity of pr and /ij, in the interior of f2. 

Lemma 6.2. Let (u k ) be the sequence of solutions of Lame system introduced in the preceding subsection, 
with N > 1. Then: 



• /il1(o,t) does not depend, in the interior of ft, on the variables Z2 and Z3; 

• /it1(o,t) does not depend, in the interior o/O, on the variable Z\. 

Corollary 6.3. Under the assumptions of lemma \6.S\ and if T is large enough: 

(80) Mi ({t€]0,T/2[})=0 

(81) M T ({te]0,T/4[, z£dQ}) =0. 

In particular, if there does not exist any boundary B -resistant ray of infinite life-length, j77| ) holds, with 
N = 1 and large enough T . 

Proof of corollary 1 6'. 51 The second part of corollary 16. 31 is an immediate consequence of (|80|l and JSTJ. 

We shall prove those two conditions by contradiction. Assume that H8U|) does not hold. Then there is, 
in the support of /x^, a point p such that: 

P ={t,z,T,(), te}0,T/2[. 

Consider D, the half-line of origin z with direction £. Let zq be the first point where D intersects the 
boundary, po — (to, z$, tq, Co) the point of the longitudinal ray coming from p whose spatial projection is 
z , and a the angle between D and the exterior normal vector to the boundary n in z (see figure 0J). If 
T is large enough (namely if ClT/2 is greater than the diameter of f2), to is in (0, T). In view of lemma 
14.51 whose assumptions are fulfilled because 1)79(1 implies (|59|l . one of the three following holds: 

• a = (corresponding to the (Hl) case of lemma |4~5*|) : 

• a = qj, = arctan(cT/ci) (corresponding to the (L — > T) case of lemma l4~5|l . 

• a = J (D is tangent to the boundary at zq, which corresponds to the (Ql) case of lemma l4~5*|l . 

Let e be an arbitrary vector, orthogonal to B. By lemma IfTTSl the following points of SZl are in the 
support of pll'- 

p e :— {t, z e , r, (), z £ :— z + ee, e > small. 
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Figure 4. Nullity of 



Consider the half-line D £ going from z £ in the direction B, and denote by zq £ its first intersection point 
with <9f2, and by a £ the angle between n £ (the exterior normal to the boundary in zq £ ), and D £ . Then: 

• if the intersection of D with dtt in p is a transverse intersection (i.e. if a 7^ then: 

Poe — J Pa, a e — J a, 

and the assumption that £1 does not have any contact of infinite order with its tangents implies 
that for small, non-zero e, a £ is close to, but distinct from a, which yields: 

(82) a e i {0,a L ,ir/2}. 

As a consequence, p £ cannot be a _B-admissible point, and D £ is not the spatial projection of a 
.B-resistant ray. For small enough e, to e is still in the interval (0,T) and thus: 

p e ^supp/xi,, 

p £ ^supp/iL (by propagation along the longitudinal flow), 
p^supp^L fbv lemma 16.21) . 

thus contradicting the definition of p; 

• if D is tangent at z to dfl, we choose e = n (which is orthogonal to B). For small e, the 
intersection of D £ and dil is a transverse intersection and is close to p. Thus, if e is small but 
strictly positive, a £ is close to, but distinct from ir/2, which again shows (|82p. and as before that 
yO is not in the support of fxi,. 

Assume now: 

#r({ie]0,T/4[, 7^0. 
Then there exists a point p such that: 

(83) p = (z,t,(,f) e supper, * e]0,T/4[, ie!l. 

Consider a transversal ray passing through p. It meets the boundary at a non-diffractive point po, at a 
certain time t$ > t, after possibly passing through diffractive points. Choose an interior point p which is 
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Figure 5. Choice of p 




Case po G Ti-T Case p € Gt 



Figure 6. Nullity of p T in the interior of 
located on this transversal ray, after all the diffractive points, but before po (see figure EJ. Thus: 

p e supppr, x en, t e]0,T/2[. 

The condition on t holds for T large enough (if suffices to take T so that the length ctT/4 of a transversal 
ray of life length T/4 is greater than the diameter of f2). 

We may choose coordinates (z 2 , z 3 ) so that ( is parallel to (0, 1,0). Let V be the plane passing through 
z and generated by the two orthogonal vectors B and £• Consider U := £1 n V, which is an open subset 
of V, and the following family of points: 

Pe ■■= (Me = z + eB,t,(), |e|<e - 

For small enough eo, z e stays in the interior of O, so that, in view of lemma 16.21 p £ is in the support 
of p*T • Let pos be the point of S£(M. x Q) where the transversal ray coming from p hits the boundary. 
Denote by n(i) the exterior unitary normal to dQ at z € 90 and, if i G "P, by n'(i) the exterior unitary 
normal to dU at i. If n(i) is not normal to the vector plane P generated by £ and B, the (non null) 
orthogonal projection of n(i) on P is parallel to n'(z). 

We shall note n , n' Q , n e , n^. instead of n(2 ), n'(zo), n(zo E ), n'(<zoe)- Two cases arise: 

• po G Wt< The point po is B-admissible, and (to being in (0,T/2)), /x^ vanishes near po by 
which implies, if \it is non null, that po is of the type Q of lemma 14.31 Consequently, no is 
orthogonal to B but not to ( (or else po should be in Gt)- Thus the (non- null) orthogonal 
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projection of no on P is ortogonal to B. Hence: 

n' ±B. 

For small enough e, zoe is close to zq. Furthermore, the preceding argument is still valid when 
replacing po by po e , so that stays orthogonal to B. Consequently dU is, in a neighbourhood 
of zq, a line segment parallel to B, contradicting the assumption that £1 has no contact of infinite 
order with its tangents. 

Po £ Gt- By the choice of p, po is not strictly diffractive and is of the Qtv*\ type of lemma l4~5l 
which implies that no is normal to the plane P. On the transversal ray coming from po, whose 
spatial projection is a geodesic curve of <9fi, n is still orthogonal to B (at least for time less than 
T). Consider points po £ , which are close to po for small e. If for one e such that \e\ < Eq, pq £ 
belongs to Tlx, we may reduce to the preceding case with p £ instead of p. Thus, we may assume: 

Ve, |e| < e =>■ p £ € Gt- 

The same argument as before shows that along rays coming from po e j the normal to the boundary 
n stays orthogonal to B. This yields a small opens subset of dil in which n is orthogonal to B, 
which shows that B is a tangent of infinite order to dtt, contradicting the assumptions on f2. 

□ 



Proof of lemma WlS. We may rewrite condition (|79|l : 



(84) VZ = 0..N, d\ +l 



-d Z2 u k - d Z3 u% 
d Zl u\ 



in iJ-^O,!? 1 ) x n). 



When N > 1, (|84|l still holds with I = 1. The characteristic manifolds of Pt and Pl are disjoint, which 
shows that (|84() holds, in the interior of fi, if one replaces u k by u\ or u\. Furthermore, on each of these 
characteristic manifolds, r does not vanish so that the operator df is elliptic. Hence: 

(85) OhH^Tjxn) 

k — > + oo 

(86) d Z2 u k L1 , d Z3 u k L1 — » in HU(0,T) x fi). 

h — ^+oc 

We have used the nullity of div ul^ to get ()85|l , and the nullity of curl u\ on the two last lines of (|84|l to 
get H86JI . We shall use properties l|85|l . (|86|l of convergence to with gain of one derivative in a classical 
way (see the proof of the propagation of defect measures), calculating the commutator of the appropriate 
operators (d zi for u\,, d Z2 and d Z3 for u\) with a "test" pseudo-differential operator of order 2. For 
example, (|85|) implies, by integration by parts: 

\fAeAl([A 7 d Zl }u k T ,u k T ) — > 

k— > + oo 

(We used that on the support of pt, v^t 1 = ||?7|| 2 .) This shows that pr does not depend upon z±. To 
prove that p^ does not depend upon Z2 nor 23, it suffices to prove the same property for the defect 
measure p^\ of v^ Ll (p is polarized along B). To do so, we use the above argument on l|8t)|) . □ 

6.3. Nullity of px on the boundary of Q,. We now assume that there exists a boundary B-resistant 
ray, 7, of infinite life length. It is transversal and its spatial projection lives on a plane curve T, contained 
in the intersection of dtt with a plane V normal to B. Furthermore, all points of 7 being gliding points, 
r is the boundary of a convex set of V. We choose T large enough so that in the interval (0,T), the 
spatial image of 7 is the entire curve T. We shall work near a point z of T, such that B is tangent at 
the order Ho = H(r) at z to dfl. This choice is made possible by the definition of Ho- Recall that by 
assumption, N > Hq. Let p be a point in the image of 7 whose spatial projection is z. 
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Notations. We choose z as the origin of the orthonormal frame (0, ei, e 2 , 63), and assume e 2 to be 
tangent, in 0, to T and e3 equal to the exterior unitary normal at to dtt. Consider the local coordinates 
(si, s 2 ) on dQ, defined by: 

Sl = Z\, s 2 = z 2 . 

The assumption on z implies that near 0: 

(87) z 3 ( Sl ,s 2 ) = s«°T (s u s 2 ), ^(si, S2 ) = Si D - 1 T 1 ( Sl , S2 ) T 3 (0,0)^0. 

The strategy of this last part of the proof is a simple one. We shall restrict condition (|79|l to the 
boundary of Q by an appropriate trace theorem. Such a theorem does not exist, in general, in spaces 
H s , s < 1/2, but in our particular case, Uj. and u\ being solutions of a differential equation which is 
transverse to the boundary, it is possible to take the trace of (|79|l on the boundary, losing as in standard 
trace theorems only half a derivative (cf proposition 17. 21 in the appendix). The boundary equations thus 
obtained will yield the nullity of ht, by a simple lemma giving bounds of L 2 norms with loss of derivatives 
flemma l6.4ll . and the usual boundary equations given by standard hyperbolic, elliptic and glancing theory. 

First step: restriction to the boundary. We shall first prove the following: 

-d Z2 u^ 2 — d Z3 uip 3 \ 

(88) curl (4 A B) ld n = B\ d Zl u k T2 WQ — > in fl£ c 1/2 ((0, T) x Oft) 

d zi u k 3 ) ^ + °° 

-d Z2 u L 2 ~ dz3 u \s , 

(89) curl («* A B) m = B 3 Zl u k L2 \ [9n —> in h£~ 1/2 ({0, T) x dQ). 



d Zl u k L3 



k- 



First of all we need to decouple condition (|79|l into one condition on the longitudinal wave and one 
condition on the transversal wave. Set: 

w k := curl (9f +1 u k A B) = io£ + w k L 

w^ L := curl(af +1 4 L AS). 

Then: 

v 2 T d 2 t w% - Aw k T = 0, v\dlw k L ~ Aw k L = 
Adding these two equations and using l|79[) we get: 

(90) 4 + v 2 L d 2 t w k L — ► in B~ 3 ((0, T) x ft). 

k — >+oo 

Furthermore, condition l|79l) twice differentiated with respect to time yields: 

(91) <9 2 ?4 + d 2 t w\ — ► in H- 3 ((0, T) x f2). 

k — >+oc 

We deduce from (|90|l and l|91|) . and i/£ being distincts: 

dfwTL — ► in H- 3 ((0,T) x fi). 

Both functions u>ji and being solutions of wave equations, proposition 17 . 21 of the appendix implies: 
curl (9f +3 4,l A B) an — > in B-J /2 ((0, T) x 90) 

Thus, if p = (t, f, 77') is a boundary point, and if f is non null, the ellipticity of 9 f iV+3 yields: 

curl (4, i AB) raOfc - T > o 0in^ 1 / 2 . 
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It remains to prove the same property when f = 0. In this case p £ £t H £l- Let M > 0. The standard 
elliptic theory (propositions 13 . 1 51 for the first line and 13.171 for the following) implies: 

(92) d Xn u% n=0 — > OinH^s 

(93) D Xn u k T{Xn=0 + E T u k TlXn=0 — > 0, in Hf 

(94) £>*»t»£ rXw= o + Siti£ riBB=0 — » OinifM, 

k— >+oc 



Sr,i G -4 1 , cri(Sr,i) = iy ||r/|| 2 - v t , l t 2 near p. 
Adding JHSJ and ijMjl. and then using (f52|. we get: 

Sr4 + Siti£ — » Oinff" 

This shows, using the Dirichlet boundary condition on u fe and the ellipticity of the operator St — 5^ 
that: 

4 >£ — > 0mH™. 

k—>+oo 

In view of l|93[l and (|94|l . the same property holds in the space H 1 on d Xn u^ L , which completes the 
proof of ® and © . 

Second step. We shall now prove that under the assumptions of subsection 16. II pyp is polarized along B: 

-KBfJ-T^B^teio.T) = /4rlte(o,T)- 
In other terms, denoting by Uj the jth component of u in the basis (ei, 63): 

/c — > + oo /c — ' + OO 

First of all, we prove the following: 

(95) a„ Uj % n fc -^ in L? a , j = 2,3. 
The sum of the last line of relations (|88|l and l|89(l yields: 

(96) d zx u\-^ 0in<- 1/2 ((0,T)x^) 



On the other hand: 

9s 1 

The Dirichlet boundary condition on u fc allows us to take out, in the preceding equality, all tangential 
derivatives. In view of l|87|) . and noting the ^--component of -J^ does not vanish at 0, we deduce from 
(I9GII that in a neighbourhood Uq of in dVt: 

(97) s «o-i an 4 0in<- 1/2 (t/ ). 

Lemma WM below allows us to take out the factor s^ 0_1 in|27| in return for the loss of Ho — 1 derivatives: 

Lemma 6.4. Let r > —1/2, p > 0, d G N* and / e H r+p (R d ), with compact support. Let y = (y 1 , ..,y<i) 
be the canonical coordinates on M. d . Then: 

||/||^<C||y?/|| H ^. 
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Proof. It is sufficient to prove the inequality with p = 1 and / € C^°(K d ). Set: 

N(f) :=-Re J (I + \ v fy Vl ^-J d v 
Cauchy-Schwarz inequality implies: 

(98) |iV(/)| < ||/|| fl -||Vl/||Hr +Ij 

A simple integration by parts yields: 



= J \f\\l + W 2 ) r - X \^ + \\ri? +rr,lj dr, 
(99) >c.||/||^, c r >0. 

To get inequality l|99|) . we used the assumption r + 1/2 > 0. Inequalities (|98|l and Q99JI yield the announced 
result. □ 

From H97fl . lemma Id^I and the assumption N — Ho > we get (|95|) with j — 2. A similar argument 
yields the same result on d n u J 3 . 

If p G 7Yl (i-e. when ct > cl) we have: 

^Iras;^ in H lan 

which implies by boundary Dirichlet condition on u k the same property on u^. 

In the elliptic case, standard elliptic theory (lemma l3.17|l yields the following boundary condition: 

Vj G {2, 3}, D n u% = Si«^ + o(l) in L?. 

This condition is still valid with T instead of L, by 1|95(1 and the Dirichlet boundary condition on u k . 

Let p! T be the defect measure of the sequence (u^ 2 , tty 3 ). Denoting by 7^3 the orthogonal projection of 
C 3 on the plane (e2, 63), one may identify p! T with 7723 /iT "23 ■ The support of p' T is, in a neighbourhood of 
p, contained in Q T and its spatial projection is contained in V. According to the propagation theorem of 
N. Burq and G. Lebeau theorem 1], there exists a function M, continuous (except possibly at hyperbolic 
points) and inversible on the support of p! T , such that M* p' T M propagates along the transversal flow 
<&t- The scalar boundary conditions above show that one may also apply the propagation theorem on 
each component an d U T3> so that M is of the form mldc2, where m is a complex valued function. 
Furthermore, /xt is polarized orthogonally to the direction of propagation, which is, along 7, the direction 
tangential to T. As a consequence, denoting by ir p the orthogonal projection of C 2 on the line generated 
by (£2, C3X the following equality yields near T: 

7T p p T = n' T , t£ (0,T) 
by scalar propagation, we get, for t € (0, T): 

Vs, 7t* t (s,p)Mt = Mr< 

This is impossible unless p! T = 0. Hence: 

4,— 0in J H|((0,T)xn), *G{2,3}. 
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Third step. We are now able to conclude to the nullity of [It- The two first lines of (|89|l may be rewritten, 
using the nullity of curl u k L : 

(d Z2 u k L1 ,d Z3 u k L1 )^0 in <- 1/2 ((0,T) x an). 
Noting that d S2 u — d Z2 u + §^-d Z3 u, we get: 

(100) d S2 u k T1 = -d S2 u k L — in H£- 1/2 ((0, T) x 90), 

k — > + oo 

which yields: 
(101) 

Consider now the first line of l|88|l . : 

(102) d zl u k Tim — » 0in< c - 1/2 ((0,T)xdO). 

We have: 

9 f a ..k dz 3 



(101) 4ir.„=o — * Oinfflg 



a n • u-JV-Ko+1/2 



ds 2 

Together with ((l"UU|) and jTQ^ this yields: 

£ (IM^) ^ ° in <- 3/2 ((0,T) x flfl). 

As a consequence, in view of the expression (I87f) of Z3, there exists a neighbourhood Uo of in dn such 
that: 

s *o-iJL ( Tl ( Sl , S2 )a„4iran) — > in <- 3/2 ((0,T) x C/ ). 
0S2 k^+00 

With lemma ROI one gets: 

(5( S i, S2 )9n4iran) - Oin <T*°- 1/2 ((0, T) x C/ ), 
which yields, using the ellipticity of the operator d S2 at p: 

Thus, since N > H : 

(103) c^^^Om^. 

Conditions (|101l) and l|103|l imply, by point b) of lemma 15. 191 the nullity of /i^i, component of /it along 
B. This complete the proof of the nullity of /j,t in view of our second step above. 

7. Appendix 

7.1. Two useful results on boundary value problems. We shall state here two lemmas concerning 
solutions of a partial diffential equation which is transverse to the boundary of an open set. The first 
one says that for such functions, the control of derivatives which are tangential to the boundary suffices 
to control all the derivatives. The second one states a trace theorem, with loss of one-half derivative, in 
all H s spaces, even if s < 1/2. As in sectional we shall work in an open subset X of R + , of the form 
X' x [0,l[, where X' is an open subset of W 1 . Let P be a diffential operator of degree r on X, of the 
following form: 

P= Qr-^, 
j=0..r 

where the Q/s are N x N matrices of tangential differential operators, with C°° [Xj coefficients, and Q 
is the identity of C^. To simplify the following statements, we suppose u S C°°(X). 
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Proposition 7.1. Let s > 0, j G N, and suppose that Pu = on X. Then: 

g C °° (X) , 3^ g C °° (X) , < C\\f>u\\ L2mHs{xl)) , 

where C does not depend on u. Likewise, if again Pu — then: 

V<peCZ°(X), 30eCg°(X), \\<pu\\l*W;h-(x>)) < C\\yu\\ H s {x) . 
Proposition 7.2. Let s € K, j E N and suppose that Pu = on X. Then: 

G C °° (X) , 3^ G C °° (X) , ||^ B «r-»=o|| lC i-x/=' < C\\0u\\ Hs{x) , 
where C does not depend on u. 

7.2. Proof of proposition 13.191 As mentionned in the introduction of section|31 whe shall assume that 
each u k is smooth enough, so that all the quantities appearing in the following calculation are well defined 
and finite. The general result may be obtained with a technical smoothing argument (cf lemma 2.8]). 
Let: 

A = A D Xn , A eA°, C k := ([P,A}u\u k ) . 

Take the support of A$ in a small enough neighbourhood of p. The operator P is formally self-adjoint. 
A simple integration by parts yields: 

C k = r (APu k , u k ) + (Au k , Pu k ) -(Au% n=0 , iD Xn u k Xn=0 ) a + (iD Xn Au k Xn=0 ,u k Xn=0 ) a , 
o 

where (., .)g is the L 2 scalar product on {x n — 0}, with respect to the measure y/g\ Xn =odx'dt. We have: 
D Xn A Q D Xn u k =A D 2 Xn u k + [D Xn ,A ]D Xn u k 

= -A Q Qu k + R Q D Xn u k 1 R Q eA° 
C k =i{A D Xn u k ,D Xn u k ) 9 - i(A Qu k ,u k ) a + (R Q D Xn u k ,u k ) a . 

So, using condition (|44|l . 

C k = i{A Q D Xn u k , D Xn u k ) d - i(AoQB^D Xri u k ,B^D Xn u k ) d - i(A Qh k , B^D Xn u k ) a 

- i{A QB_ 1 D Xn u k 1 h k ) d + (R D Xn u k , B- X D Xn u k ) a + o(l), k -> +oo 

Take Ao of the form Tq To, where To G Ao is scalar, elliptic at p and has support in a small neighbourhood 
of p. Then: 

C k = i(E D Xn u k , D Xn u k ) d - i(T Qh k , B^T D Xn u k ) - i(QB^T D Xn u k ,T h k ) + o(l), k -> +oo 
Eq g A , Eq = TqTq — _B1 1 TqToQS_i + -B^jiJo- 
Denoting by to the principal symbol (which is scalar) of To, we have: 

cto(£o) = |*o | 2 (1- tUb-ita) ■ 
Since qi vanishes at p, we may choose to with support in a small enough neighbourhood of p such that: 

<j(Eo) > l/2\t Q \ 2 

(in the sens of quadratic positive hermitian forms). The weak Garding inequality, applied to the operator 
E Q - 1/2T *T , thus yields: 

liminf Rje(E D Xn u k ,D Xn u k ) - l\\T D Xn u k \\ 2 L2 > 0. 

k^ + oo 2 n a 

This implies, using the convergence to of h k in H~: 

liminf ImC fe > -\\T D X u k \\ 2 T 2. 

k^ + oo 4 a 
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Thus D Xn u^ Xn=0 is bounded in L 2 ~ d , which yields, with the boundary condition, that ut.j is bounded 
in H~ g . The proof of is complete. 
When /i is null near p, we have: 

hm C ={n, )=0, 

which yields iJSDJ. 

Point b) of proposition 13.191 may be seen as a consequence of the propagation theorem of [2] . The 
assumptions imply that any uniform Lopatinsky boundary condition holds on the traces of u k , which 
shows that the measure /i propagates near p with any smooth multiplicative factor, which is impossible 
unless /i is null. 
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